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Regional  Analysis  of  Low  Flow  Characteristics 
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(e.g.  Technical  Transfer  Conference) . 

Overall,  the  results  of  the  various  methods  for  regionalizing  low  flows 
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investigations  in  the  literature.  However,  this  work  has  identified 
two  main  areas  of  additional  investigation  which  we  believe  should  be 
pursued: 

a)  Detailed  Assessment  of  Trend,  Persistence  and  Randomness 

b)  Extension  of  Regional ization  to  Other  Regions 

Thank  you  for  the  opportunity  to  have  undertaken  this  most  interesting 
investigation  on  your  behalf. 

Yours  very  truly 

CUMMING  COCKBURN  LIMITED 


H.  S.  Belore,  P.  Eng. 

Manager,  Water  Resources  Division 
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ABSTRACT 

REGIONAL  ANALYSIS  OF  LOW  FLOW  CHARACTERISTICS 


The  identification  of  suitable  low  flow  characteristics  within  a  water- 
course is  most  easily  accomplished  using  continuous  hydrometric  data 
recorded  for  the  stream.  However,  this  is  limited  to  the  availability 
of  suitable  long  term  discharge  records  at  the  location  of  interest. 
The  objective  of  this  investigation  is  to  develop  a  practically  oriented 
technique  for  estimating  low  flow  characteristics  at  ungauged  locations. 

Initially,  a  background  review  was  undertaken  to  confirm  the  most  appro- 
priate statistical  techniques.  It  was  confirmed  that  the  Gumbel  III 
extreme  value  distribution  is  acceptable  for  analysis  of  single  station 
low  flow  data.  However,  as  an  alternative  distribution  for  samples  with 
large  negative  skewness,  it  was  determined  that  the  Three  Parameter  Log 
Normal  Distribution  should  be  used. 

The  background  review  also  confirmed  that  several  previous  investigators 
utilized  multiple  linear  regression  techniques  to  regionalize  low  flow 
characteristics.  On  the  basis  of  the  literature  review,  it  was  decided 
to  -include  the  following  parameters  in  the  development  of  a  physiograph- 
ic data  base  for  the  Southwestern  and  West  Central  Regions  of  Southern 
Ontario:  drainage  area  (DA),  mean  annual  snowfall  (MAS),  mean  annual 
precipitation  (MAP),  mean  annual  runoff  (MAR),  area  controlled  by  lakes 
and  swamps  (ACLS),  Latitude  (LAT),  Longitude  (LONG),  10/85  slope  (SLP), 
stream  length  (LEN),  years  of  daily  data  (DYRS),  and  Base  Flow  Index 
(BFI).  After  preliminary  analysis  and  screening,  a  data  base  was 
developed  for  a  total  of  65  watersheds. 

Several  alternative  methods  were  considered  for  regionalizing  low  flow 
characteristics  (i.e.  7020  and  O95)  including;  i)  Regressional  Analysis; 
ii)  Index  Method;  iii)  mapping  proration  and  iv)  mapping  isolines.  The 
methods  developed  were  then  applied  and  evaluated  for  10  test  stations 
to  confirm  the  usefulness  in  predicting  low  flow  characteristics  for 
ungauged  watersheds. 

Regressions  utilizing  drainage  area  (DA),  base  flow  index  (BFI)  and 
stream  length  (LEN)  as  independent  parameters  were  found  to  provide 
reasonable  estimates  of  low  flow  values.  The  index,  or  graphical  method 
was  also  found  to  give  reasonable  estimates  of  low  flows.  Both  these 
methods  provide  low  flow  estimates  superior  to  those  obtained  solely 
from  area  proration  to  gauged  sites. 

Some  requirements  for  additional  investigations  were  also  identified, 
including  an  examination  of  possible  reasons  for  persistence,  trend  and 
non-randomness  in  the  available  low  flow  data  set  and  further  extension 
of  the  index,  mapping  and  regression  estimating  methods  to  other  regions 
in  the  Province. 
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1.0  INTRODUCTION 

1.1  General 

The  knowledge  of  the  hydrologic  conditions  which  exist  during  low 
flow  conditions  can  be  of  primary  importance  for  some  water- 
courses. For  example,  when  analysing  water  quality  conditions, 
the  low  flow  conditions  within  the  watercourse  are  a  major  concern 
to  both  the  landowner  and  the  Ontario  Ministry  of  the  Environ- 
ment. Other  uses  of  low  flow  information  may  include  the  follow- 
ing: 

i)  In-stream  pollutant  analyses  (point  and  non-point  sources) 
ii)  Reservoir  design  (low  flow  augmentation) 
iii)  Environmental  appraisals 
iv)  Feasibility  of  small  hydro  developments 
v)  Water  supply  and  evaluation  for  water  taking  permits 
vi)  Base  flow/groundwater  recharge  and/or  contamination 
analysis 
vii)  Stream  fisheries  assessments 
viii)  Analyse  effects  of  changes  in  watershed  on  low  flows 
(e.g.  deforestation,  urbanization) 
ix)  Agricultural 

The  identification  of  low  flow  characteristics  within  a  water- 
course is  most  easily  accomplished  using  continuous  hydrometric 
data  recorded  for  the  stream.  However,  the  effectiveness  of  this 
method  is  limited  since  a  hydrometric  recording  gauge  may  not  be 
located  in  the  vicinity  of  the  study  site. 

Another  available  source  of  information  pertaining  to  drought  con- 
ditions is  the  "Low  Flow  Characteristics"  maps  which  were  publish- 
ed several  years  ago  (more  than  10  years  in  most  cases)  by  the 
Ministry  of  the  Environment.  These  maps  were  produced  for  most 
regions  of  Ontario  and  include  exceedance  and  extreme  values 
describing  low  flow  characteristics  of  gauged  watersheds. 
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However,  a  significant  amount  of  additional  streamflow  measure- 
ments are  now  available  across  the  Province.  Furthermore,  diffi- 
culties can  arise  in  transferring  this  information  to  ungauged 
sites,  since  generalized  procedures  are  not  available. 

Therefore,  this  project  was  conceived  to  develop  the  methodology 
for  a  regional  estimating  technique  for  estimating  low  flow  char- 
acteristics at  ungauged  locations.  The  test  areas  selected  were 
the  Southwestern  and  West  Central  regions  of  the  Province  of 
Ontario  (see  Figure  1.1). 

Regional  hydrologic  analyses  can  take  several  forms  and  can  util- 
ize regression  equations,  index  methods,  mapped  isolines,  etc. 
However,  all  of  these  techniques  serve  to  correlate  selected  char- 
acteristics of  an  existing  hydrometric  data  base  to  physiographic 
and  meteorologic  parameters  which  can  be  measured  at  both  gauged 
and  ungauged  watersheds.  By  utilizing  the  relationship  developed 
using  the  gauged  watersheds,  the  required  characteristics  of 
streamflow  can  be  estimated  for  any  ungauged  watershed  within  the 
region. 

1.2    Study  Objectives 

The  main  objectives  are  to  develop  a  technique  for  providing  esti- 
mates of  low  flow  characteristics  for  a  watercourse  based  on  the 
physical  parameters  of  the  watershed  and  appropriate  meteorologic- 
al variables.  This  includes  a  demonstration  application  of  the 
technique  within  the  region  of  development. 

The  technique  could  then  be  further  developed  and  adapted  in  order 
to  provide  estimates  of  low  flows  for  ungauged  watersheds  at  other 
locations  in  the  Province. 
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2.0  LITERATURE  REVIEW 

2.1  General 

A  summary  of  the  relevant  results  of  the  review  of  the  various 
investigations  is  given  in  Table  2.1. 

The  following  sections  summarize  highlights  from  selected  refer- 
ences. The  reader  should  refer  to  the  original  reference  before 
attempting  to  apply  specific  equations  since  the  units  and 
parameter  definitions  vary  from  report  to  report.  The  intent  here 
IS  only  to  identify  important  parameters  for  further  consideration 
in  this  investigation. 

Median  Drought  Flows  at  Unqai.gPri  Sites  in  .Sonthprn  n......   (pi^,,^ 

P.  and  Condie,  R.,  1986) 

A  methodology  was  suggested  and  explored  for  predicting  median 
drought  flows  at  ungauged  sites  in  Southern  Ontario.  Minimum  aver- 
age flows  for  durations  of  3,  10,  14,  21  and  28  consecutive  days 
were  determined  each  year  in  the  time  period  June  1  to  May  31  for 
74  stations  in  Southern  Ontario.  A  regression  equation  was 
developed  for  predicting  the  median  10-day  low  flow  using  drainage 
area  (DA),  mean  annual  precipitation  (MAP),  and  base  flow  index 
(BFI).  The  precipitation  parameter  was  not  significant  and, 
therefore,  not  included  in  the  equation.  The  standard  error  of  the 
low  flows  was  found  to  be  37%  and  the  equation  was  in  the  form: 

MED  (10)   =  0.000075  DA  +  1.515  BFI  -  0.356 

The  median  10-day  low  flow  was  then  used  to  produce  a  slope 
relationship  to  relate  other  low  flow  durations  to  the  MED(IO) 
value.  This  equation  is  expressed  as: 

GRADMEN  =  2.632  MED(IO)  +  0.46399 
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where  GRADMEN  is  the  slope  of  the  median  flow-duration 
relationship. 

A  problem  noted  by  the  authors  is  the  determination  of  BFI  for 
ungauged  sites.  The  parameter  is  used  as  an  index  of  catchment 
geology. 

Modelling  of  Low  Flow  Frequency  Distributions  and  Parameter  Estima- 
tion  (Condie,  R.  and  G.  Nix.  1974) 

The  fit  of  four  theoretical  frequency  distributions  was  examined  in 
this  study.  The  Gumbel  III  distribution  was  selected  as  the  best 
fit  for  extreme  value  analysis  of  low  flows  since  it  produced  33 
acceptable  fits  out  of  a  possible  38  for  the  Canadian  rivers  used 
in  this  study.  A  program  FLOINT  was  written  to  handle  low  flow 
values  of  zero  using  a  joint/conditional  probability  function.  A 
discussion  of  the  lower  boundary  limit  being  between  zero  and  the 
lowest  recorded  flow  concluded  that  the  lower  boundary  can  be  less 
than  zero  since  it  represents  the  curve  of  the  line  and  not  the 
actual  flow. 

Low  Flow  Frequency  Analysis  (Program  LOFLOW)  (R.  Condie  and  L. 
Cheng,  1983) 

The  Gumbel  Type  III  distribution  is  discussed  as  well  as  conceptual 
problems  found  from  the  Condie,  Nix  1975  study.  It  was  found  that 
The  Gumbel  III  distribution  becomes  unstable  if  the  sample  has  a 
skewness  of  less  than  -1.08.  This  study  concluded  that  the  Gumbel 
III  distribution  should  be  used  as  the  basic  method  for  low  flow 
frequency  analysis.  As  an  alternative  for  samples  with  large 
negative  skewness,  it  was  suggested  that  a  3-parameter  log  normal 
distribution  should  be  used. 

The  only  criteria  of  successful  fit  to  the  Gumbel  III  is  that  the 
lower  boundary  parameter  should  not  be  larger  than  the  smallest 
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member  of  the  sample.  Parameter  estimation  should  proceed  in  the 
order  of;  maximum  likelihood,  smallest  observed  drought  and  by 
moments.  The  program  LOFLOW  was  written  to  analyse  samples  of  low 
flows  taking  these  concerns  into  account. 

Computation  of  the  Base  Flow  Index  Proor.m  (RPT.nrv^  (p.  5^3^  and 
R.  Condie,  1983) 

It  was  determined  that  low  flows  are  very  dependent  on  catchment 

geology,  and  the  Base  Flow  Index  was  suggested  as  a  practical 

technique  to  provide  an  index  of  this  parameter.  A  program  was 

written  to  calculate  a  base  flow  index  (BFI)  as: 

BFI  =  r-unoff  under  base  flow  separation  linp  yr. 
total  runoff  yp. 

The  required  data  was  obtained  from  records  published  by  the  Water 
Survey  of  Canada,  and  averaged  BFI  values  were  computed  over  a 
period  of  several  years. 

The  main  problem  with  this  technique  was  determined  to  be  the 
requirement  for  discharge  measurements.  However,  it  was  pointed 
out  that  the  BFI  parameter  was  found  to  be  fairly  stable  over  botn 
wet  and  dry  years  and,  it  was,  therefore,  concluded  that  it  could 
be  estimated  from  as  little  as  one  to  two  years  of  gauged  record. 

Estimates  of  low  Flows  Using  Watershed  Climatic  ParametPr.  (Boyer, 
P.G.  and  M.  Chang,  1977) 

A  regional  analysis  was  undertaken  for  12  unregulated  mountainous 
watersheds  with  drainage  areas  ranging  from  164  km2  to  2345  km2 
(64-916  mi2  located  in  West  Virginia).  There  was  considerable 
effort  put  forward  in  developing  physical  watershed  parameters  for 
which  the  evaporative  power  and  drought  potential  were  related 
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logically.  Eighteen  watershed  parameters  were  tested;  drainage 
area,  length  of  perimeter,  mean  elevation,  mean  latitude,  maximum 
relief,  relief  ratio,  basin  slope,  form  factor,  main  channel  slope, 
main  channel  length,  drainage  density,  stream  frequency,  stream 
order,  average  angle  of  tributaries,  aspect,  length-width  ratio, 
coefficient  of  compactness  and  percent  of  forest  cover.  Climatic 
parameters  used  for  analysis  were  7-day  mean  maximum  temperature 
and  maximum  consecutive  days  without  precipitation  for  seasonal 
periods  of  August  to  October. 

The  final  model  used  3  watershed  and  2  climatic  parameters,  namely; 
watershed  perimeter  (CO),  the  main  channel  length  in  miles  (LEN), 
the  Watershed  Form  Factor  (WF)  (area/maximum  watershed  length^). 
The  7-day  10  year  maximum  temperature  (Tiq),  a  function  of  mean 
elevation  and  latitude,  and  the  September  10  year  maximum  consecu- 
tive rainless  days  (Rio),  and  finally,  a  function  of  watershed  mean 
elevation.  A  prediction  equation  was  developed  for  the  1:10  year 
7-day  low  flow  with  a  standard  error  of  estimate  of  about  30%  of 
the  observed  mean: 

In  (70io)  =  37.0469  +  0.03886(CO)  -  0.0401(LEN)  -  3.6743(WF) 
-  .5771  Tio  +  1.1665  Riq 

It  was  concluded  that  use  of  climatic  parameters  which  index 
drought  conditions  can  improve  low  flow  estimates,  and  that  mean- 
ingful low-flow  estimates  can  be  obtained  from  climatic  and  water- 
shed parameters  or  watershed  paramaters  alone  in  a  mountainous 
humid  region. 

Analysis  of  Low  Stream  Flows  on  Cape  Breton  Island,  Nova  Scotia 
(Environment  Canada,  1978) 

This  investigation  included  frequency  analysis,  using  Gumbel  III, 
for  1,  3,  7,  10,  30  day  low  flow  durations  based  on  11  hydrometric 
stations  with  drainage  areas  ranging  from  216  km^  to  3.6  km^  (142 
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mi2  -  2.4  ml2).  Stations  with  periods  of  record  from  10  to  59 
years  were  used  and  those  with  man-made  regulation  were  deleted. 

It  was  determined  that  extreme  low  flow  estimates  for  the 
short-term  hydrometric  stations  may  be  high.  Seasonal  analysis 
found  higher  low  flows  for  the  sunmer  period.  It  was  suggested 
that  this  was  due  to  the  climatic  and  physical  location  of  the 
stations,  and  it  was  noted  that  evaporative  power  may  also  be 
important.  As  well,  a  physiological  data  base  was  prepared  for  a 
future  regression  study  of  13  parameters  relating  to  the  stations 
examined  (see  Table  2.1). 

Journal  of  the  Hydraulic  Division  -  Characteristics  of  Low  Flows 
(Task  Committee,  Surface  Water  Hydrology,  1980) 

This  paper  summarizes  the  findings  of  a  U.S.  committee  regarding 
the  types  of  low  flow  information  needed,  various  low  flow  charac- 
teristics and  their  accuracy,  and  suggests  further  analysis  and 
data  collection  needs.  The  main  topics  were  discussed  according  to 
gauged  sites  and  ungauged  sites. 

Gauged  Site 

The  most  widely-used  index  of  low  flow  in  the  U.S.  is  the  7-day 
10-year  low  flow  (70io)-  However,  it  was  also  found  that  other 
durations  of  from  1-273  days  have  been  based  on  data  records  of 
15-20  years  or  more. 

Ungauged  Sites 

It  was  found  that  low  flow  characteristics  depend  largely  on  the 
geology  of  the  basin  and  on  losses  from  evapotranspiration.  For 
example,  the  incorporation  of  basin  rock  type  into  regional  analy- 
sis of  low  flows  increased  prediction  accuracy  for  equations  devel- 
oped in  Virginia.  A  regional  low  flow  forecasting  model  developed 
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(by  Wright,  1970)  for  southeast  England  relates  lowest  mean  flow 
for  the  forecast  year  to,  drainage  area,  mean  annual  runoff,  geol- 
ogy index,  mean  summer  catchment  rainfall  in  the  same  year  as  fore- 
cast and,  mean  winter  catchment  rainfall  (October-March)  preceding 
forecast  low  flow. 

Effects  of  Man  on  Low  Flows  (Riggs,  H.C.) 

This  investigation  examined  the  climatic,  geologic  and  topographic 
influences  on  low  flow  characteristics  in  view  of  man's  influence. 
Man-made  changes  that  affect  low  flows  were  identified  as  changes 
in  vegetation,  drainage  practices,  diversions,  return  flows, 
groundwater  pumping,  sewage  treatment  plants  and  other  effluent 
di  scharge. 

It  was  pointed  out  that  the  aquifer  characteristics  are  also 
extremely  important  to  low  flows.  It  was  postulated  that  evapo- 
transpi ration  and  changes  to  ground  cover  can  likely  be  related  to 
low  flow  increases.  Urbanization  is  not  a  significant  factor  since 
pluses  (i.e  sewer  leakage)  and  minuses  (i.e.  reduced  infiltration) 
tended  to  even  out.  Reservoir  regulation  usually  increases  low 
flows  due  to  minimum  outflow  requirements.  Diversions  are  highly 
basin  dependent  and  are  difficult  to  assess. 

It  was  noted  that  swamps  apparently  had  very  small  low  flows  in 
extreme  drought  cases.  Therefore,  it  was  postulated  that  changes 
to  swamps  (e.g.  by  infilling,  etc.)  might  not  have  a  significant 
impact  on  downstream  low  flow  characteristics.  Pumping  can  have  a 
large  effect  depending  on  duration  location  and  basin  geology.  The 
three  main  parameters  for  low  flow  are  related  by;  climate  provid- 
ing input  and  output,  geology  and  topography  influencing  infiltra- 
tion to  aquifers,  and  then  geologic  parameters  affecting  transport 
to  streams. 
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Man's  influence  on  precipitation,  temperature  and  irrigation  are 
not  likely  to  effect  low  flows  while  changing  vegetation  cover 
(evapotranspiration) ,  regulation  and  diversion  are  likely  to  have 
the  most  evident  effects  on  low  flow. 

Hydrology  of  7-Dav  10  Year  Low  Flows  (Singh,  K.P.  and  J.B.  Stahl , 
1974) 

This  paper  discusses  three  pertinent  factors  in  the  broad  headings 
of;  1)  basin  factors;  2)  climatic  factors,  and  3)  sociological 
factors.  The  main  discussion  is  based  upon  man's  influence  and 
requirement  for  low  flow.  It  is  suggested  that  basin  input  to  a 
7Qio  low  flow  study  includes  information  on  daily  stream  flows, 
soil  and  stream  characteristics,  physiographic  and  hydrologic 
diversion  boundaries,  municipal  and  industrial  water  use  and/or 
man-made  lakes  and  flow  regulation.  Again,  it  is  stressed  in  the 
conclusions  that  the  magnitude  of  the  drainage  area  is  the  most 
important  factor  affecting  the  natural  70io  of  hydrological ly 
simi lar  basins. 

An  Infiltration  Index  Useful  in  Estimating  Low  Flow  Characteristics 
of  Drainage  Basins  (Armbruster,  J.T.,  1976a) 

This  investigation  concluded  that  the  basin  geology  was  an  over- 
riding factor  in  describing  low  flow  characteristics.  Using  an  SCS 
soil  type  based  infiltration  index,  the  author  of  this  study  was 
able  to  improve  standard  errors  in  existing  prediction  equations. 
The  soil  index  was  determined  by  percentage  of  basin  of  soil  type 
A,  B,  C  or  D  and  applying  a  weighting  of  10,  5,  1,  1  respectively 
and  summing  for  the  entire  basin.  The  drainage  area  and  a  precipi- 
tation index  were  also  found  to  be  significant  factors  in  estima- 
ting low  flows.  The  final  equation  takes  the  following  form: 

log  70io  =  -4.658  +  1.068  log  DA  +  2.183  log  MAP  +  .218  SI 
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Technical  Manual  for  Estimating  Low  Flow  Frequency  Characteristics 
of  Streams  in  the  Susquehanna  River  Basin  (Armbruster,  J.T.,  1976b) 

This  technical  report  adopted  the  infiltration  index  described  in 
the  earlier  reference  by  Armbruster  which  examined  step-by-step 
procedures  for  low  flow  frequency  characteristics.  In  addition, 
distinctions  were  made  between  large  and  small  drainage  basins. 
Large  basins  have  drainage  areas  of  larger  than  5,120  km^  (2,000 
mi2).  The  parameters  used  in  the  regressions  were;  drainage  area, 
mean  annual  precipitation,  basin  surface  storage,  index  of  relative 
infiltration.  As  well,  an  index  coefficient  of  the  1-day  flow  to 
the  3-day  flow  of  .93  was  determined.  For  large  basins,  a  graph  of 
correction  factors  for  channel  length  was  produced  and  applied  to 
the  regression  equation  as; 

log  (7020)  =  -2.029  +  1.052  log  DA  +  .753  log  (MAP)  +  CF 

where  CF  is  the  correction  factor  from  the  graph. 

Combining  Estimates  of  Low  Flow  Characteristics  of  Streams  in 
Massachusetts  and  Rhode  Island  (Tasker,  G.D.,  1975) 

This  report  presents  a  statistical  analysis  of  baseflow  measure- 
ments and  a  discussion  of  the  accuracy  produced  by  varying  the 
sample  size.  It  was  concluded  that  empirically  estimated  values  of 
variance  of  the  estimated  70io  low  flow  from  its  true  value  indi- 
cate that  more  than  six  or  eight  baseflow  measurements  add  little 
to  the  confidence  with  which  such  an  estimate  is  made.  As  well, 
the  number  of  baseflow  measurements  needed  to  meet  a  specific 
accuracy  of  low  flow  characteristics  can  be  significantly  reduced 
by  considering  prior  knowledge  of  the  low  flow  characteristics 
quantitatively.  The  formula  is: 

q  ^  qr  +  (Sr2/Sb2)  qfa 
1  +  (Sr2/Sb2) 
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where   q  is  logarithm  of  estimated  low  flow  characteristics 
Qf-  is  log  Qp  flow  estimated  by  regression  of  basin 

parameters 
Qb  is  log  Ob  flow  estimated  by  regression  of  base  flow 
measurements 
Sb2  and  Sr^  are  variances  in  square  log  units  of  estimates  from 
true  value  of  low  flow  characteristics 

Probability  Distribution  of  Low  Flows  (Matalas,  N.C.,  1963) 

This  study  examined  the  fit  of  four  probability  distributions  using 
both  the  method  of  maximum  likelihood  and  moments.  The  Gumbel  and 
Pearson  Type  III  distributions  fitted  the  data  equally  well  and 
were  more  representative  of  the  probability  distribution  than  the  3 
Parameter  Log  Normal  or  Pearson  Type  V  distributions.  Goodness  of 
fit  was  based  on  observed  minimum  flow,  the  lower  limits  of  the 
four  probability  distributions  and  the  relation  between  skewness 
and  kurtosis  of  the  low  flow  data.  It  was  recommended  that  future 
statistical  studies  of  low  flows  should  use  the  method  of  maximum 
likelihood  estimates  of  parameters  rather  than  moment  estimates. 

Drought  Flows  and  Receiving  Water  Assessment  (Logan,  L.,  1986) 

This  study  examined  low  flows  in  an  environmental  context  with 
respect  to  waste  water  and  failures  to  satisfy  water  quality 
requirements.  Low  flow  design  estimates  for  Ontario  Central  - 
Southwest  were  reported  to  be: 


Central      Southwest 
(m^/s/km^)     (m3/s/km^) 


7020  1-0    -     0.8    ) 

70io  1.1    -     0.9    )   ±.3  m3/s/km2 

702  1-8     -     1.3     ) 
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The  M.O.E.  guideline  on  drought  flow  design  is  7020«  This  implies 
that  lower  low  flows  will  cause  water  quality  violation  in  BOD  and 
other  wastes. 

Flood  Prone  Area  and  Low  Flow  Analysis  Study  in  the  Kawartha  Region 
(MacLaren  Plansearch,  1981) 

This  study  developed  a  regression  equation  for  predicting  low  flow 
estimates  in  the  Kawartha  region.  The  regression  equations  took 
the  form: 

§^  =  1.64  X  10-3  SI3.171  ALS-0.116 
g|  =  6.90  X  10-4  sil.365  aLS-0-102 


where  Qi  =  average  minimum  daily  flow 

Qs  =  average  minimum  summer  flow 
DA  =  drainage  area 

SI  =  index  of  average  soil  moisture  holding  capacity 
ALS  =  %   area  covered  by  swamp  or  wet  land 

Approach  for  Frequency  Analysis  of  Multiyear  Drought  Durations 
(Lee,  K.S.,  Sadeghipour,  J.,  and  J. A.  Dracup,  1986) 

An  approach  for  frequency  analysis  of  multiyear  drought  durations 
was  presented.  The  problem  of  sample  size  is  discussed  with  a 
technique  for  smoothing  frequency-curve  irregularities  of  drought 
durations.  Piecewise  linear  regressions  for  the  linearized  termin- 
ation probability  curve  define  the  parameters  of  the  logistic 
model.  The  model  then  can  be  used  to  determine  the  return  periods 
of  multilinear  droughts.  This  report,  although  interesting,  is  not 
directly  applicable  to  the  present  investigations. 
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Determining  Streamflows  from  Geomorphic  Parameters  (Osborn,  J.F., 
1974) 

Osborn  examined  the  many  interconnected  relationships  of  geomorphic 
parameters  of  basins  to  flows.  An  example  of  the  basin  character- 
istics interrelationship  is  that  LS  =  1.4  (DA)0'6,  where  LS  is 
stream  length  and  DA  is  area.  It  is  noted  that  the  coefficient  and 
power  values  vary  for  different  geographic  areas. 

The  equations  for  estimating  average  annual  flow  (QAA)  took  the 
form  QAA  =  C"  (P)^  (DA)^  where  P  is  average  annual  precipitation 
and  the  coefficient  and  powers  are  determined  from  regional 
regressions. 

It  is  noted  that  minimum  or  low  flows  tend  to  be  more  sensitive  to 
a  larger  number  of  parameters  than  flood  flows.  Assuming  linear 
flow-frequency  relationships,  a  parameter  07L1P  was  developed  to 
relate  the  relative  maximum  size  of  low  flow  in  one  stream  to 
another.  This  value  was  then  regressed  with  basin  parameters  to 
establ ish  702. 


^^2  -  e  L      300{P)   -■ 


where  B  is  a  coefficient 
H  is  basin  rel ief 
P  is  slope  of  recurrence  interval  graph 

The  P  parameter  also  represents  a  measure  of  the  stability  of  the 
natural  basin  storage  from  which  the  low  flow  is  derived.  The 
value  702  ^^s  then  used  as  an  index  flow  to  predict  other  recurr- 
ence interval  flows  as: 


7Q20  =  !ii!l 


Where  P  is  slope  and  C  is  a  coefficient.  For  example, 
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log  7020  =  log  (P(7Q2))  -  log  C 

Relationships  were  also  found  between  peak  flows  and  low  flows.  It 
was  determined  that  if  basins  were  of  the  same  average  annual  run- 
off then  you  assume  C  is  essentially  the  same  and  examine  the  peak 
flows,  if  they  are  different,  then  low  flows  must  be  different  and 
the  basins  are  not  hydrological ly  similar. 

Hydrologic  Design  Methodologies  for  Small  Scale  Hydro  at  Ungauged 
Sites  (Acres,  1985) 

This  study  developed  a  regression  relationship  between  basin  physi- 
ographic parameters  and  "turbinable  flow"  which  was  related  to 
exceedance  curves.  Exceedance  curves  were  non-dimensional ized  by 
dividing  flow  values  by  the  mean  annual  flow.  For  ungauged  sites, 
a  regression  equation  provided  estimates  of  mean  annual  flow  from 
values  of  MAR  and  area. 

MAF  =  0.0000369  MAR-DA 

This  was  then  related  to  a  regional  indexed  curve  to  determine 
values  of  exceedance  by  examining  the  regional  mean  annual  flow  to 
the  sight  mean  annual  flow. 

As  well,  two  other  methods  were  compared,  proration  and  ARMA 
(Autoregressive  Moving  Average  Methodology).  It  was  determined 
that  proration  provided  the  best  practical  method  for  determining 
discharge.  The  proration  formula  took  the  form: 

Oijk  =  QIijk(DA  MAR)   /  (DAI  MARI) 

where  i  =  year,  j  =  month,  k  =  day 

The  ARMA  method  was  used  to  synthesize  daily  and  monthly  flows  from 
required  parameters  within  the  basin  over  a  time  series. 
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Low  Flow  Studies  (Institute  of  Hydrology,  1980) 

The  Institute  of  Hydrology  (U.K.)  analysed  low  flows  using 
regression  techniques.  An  innovative  feature  was  the  development 
of  a  geologic  index  called  the  baseflow  index.  This  index  can  be 
determined  from  streamflow  measurements  (for  equation  see  Condie, 
1986).  In  ungauged  areas,  comparative  area!  methods  can  be  used  to 
estimate  the  BFI  or  it  can  be  estimated  from  underlying  geology. 

Low  flow  duration  curves  were  regressed  from  basin  parameters. 
Then  a  comparative  analysis  of  the  change  of  slope  for  different 
d-day  durations  allowed  for  a  prediction  of  any  d-day  low  flow.  In 
the  study,  much  work  was  done  to  provide  non-dimensional  graphs  and 
to  linearize  them  for  simplifying  the  estimation  procedures. 

The  regional  equations  take  the  form: 

Q95  (10)  =  Ci  BFI  +  C2  MAP  -  C3 
and  Q95  (10)  =  C4  BFI  +  C5  LEN  -  Cs 

where  MAP  =  mean  annual  precipitation 
LEN  =  stream  length 

for  region  five. 

As  well,  the  slope  is  regressed  from  parameters  as: 

logio  GRADMEN  Q95  =  C7  MAP  -  Cs  Q95(10)  -  Cg 

These  two  regression  equations  were  then  used  to  estimate  095  (d) 
for  any  d-day  duration  by  the  formula: 

O95  (d)  =  Q95  (10)  [1  +  (D-10)  GRADMEN  O95]. 

Similar  relationships  were  developed  for  mean  annual  minimum  10-day 
low  flow  MAM  (10).  Included  in  the  report  describing  this  five 
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year  study  is  a  fairly  extensive  literature  review  which  is  para- 
phrased below. 

Review  of  the  Institute  of  Hydrology  Background  Review  (Institute 
of  Hydrology,  1980) 

"Many  examples  worldwide  employ  catchment  characteristics  indexing 
catchment  size,  shape  and  climate  and  most  refer  to  the  additional 
need  to  consider  the  geology  of  the  catchment  in  order  to  explain 
more  of  the  variability  of  low  flows.  Wright  (1970)  presents  a 
prediction  equation  for  the  lowest  flow  in  a  year  using  catchment 
slope  and  area  as  Independent  variables.  Riggs  (1973),  Wesche  et 
al  (1973),  and  Skelton  (1974),  make  use  of  the  cross-sectional 
properties  of  streams  to  correlate  low  flows  with  neighbouring 
gauged  catchment  data.  For  the  United  Kingdom,  general  information 
on  behaviour  of  British  rivers  is  given  by  (Ward,  1968)  and  also 
(Rodda,  1976).  The  interaction  between  surface  and  groundwater  for 
major  aquifers  is  discussed  in  Ineson  et  al  (1965).  Much  design 
work  has  been  based  on  the  lowest  recorded  flow.  Flow  duration 
curves  and  their  applications  in  the  U.K.  are  given  in  Hoyle 
(1963).  They  are  discussed  by  Boulton  (1965)  where  the  concept  of 
a  Minimum  Acceptable  Flow  is  described  and  the  factors  likely  to 
enter  into  its  calculations  are  considered." 

Regional  Flood  Frequency  Analysis  for  Ontario  Streams,  Volume  2, 
Multiple  Regression  Method  (Moin,  S.M.A.,  and  M.A.  Shaw,  1986 

This  report  includes  a  discussion  of  the  Base  Flow  Index.  The  BFI 
was  calculated  and  mapped  for  2  flood  regions  identified  in  the 
Province. 

The  regression  equations  with  BFI  for  predicting  flood  flows  pro- 
vided predictions  superior  to  those  derived  using  the  secondary 
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regression  equations  without  the  BFI  parameter.  The  latter  was 
found  to  be  second  only  to  drainage  area  as  a  significant  para- 
meter in  describing  variance  of  flood  flows. 

The  BFI  was  found  to  be  extremely  sensitive  to  basin  storage  condi- 
tions and,  therefore,  could  not  be  used  in  the  primary  regression 
equations  for  the  northern  region  due  to  the  large  number  of  lakes 
and  swamps. 

Correlations  with  SCS  soil  groups  indicated  that  BFI  could  be 
approximated  as  follows: 

Soil  Group  A  >  0.60± 

Soil  Group  B  0.45± 

Soil  Group  C  0.30± 

Soil  Group  D  <  0.25± 

2.2    Summary  of  Literature  Review 

In  summary,  it  was  found  that  many  investigations  have  produced  low 
flow  durations  from  1  to  273  days  in  length,  for  return  periods  of 
1  to  100  years.  In  general,  it  has  been  determined  that  the  Gumbel 
III  distribution  best  fits  low  flow  data  for  gauged  streams 
(Condie,  Nix,  1975,  Matalas,  1965  and  Table  2.1).  On  the  other 
hand,  very  few  investigations  have  considered  techniques  for  devel- 
oping regionalized  flow  duration  curves. 

Parameters  found  to  be  significant  in  regression  equations  general- 
ly tended  to  be  drainage  area,  base  flow  index,  and  mean  annual 
precipitation.  However,  some  parameters  that  have  not  been  highly 
correlated  in  one  region  were  found  to  be  important  in  other 
regions  (Institute  of  Hydrology,  1980). 

From  the  available  literature,  it  was  decided  to  include  the  foll- 
owing parameters  in  development  of  the  physiographic  meteorologic 
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and  hydrometric  data  base;  drainage  area,  (DA),  mean  annual 
snowfall  (MAS),  mean  annual  precipitation  (MAP),  mean  annual  runoff 
(MAR),  area  controlled  by  lakes  and  swamps  (ACLS),  Latitude  (LAT), 
Longitude  (LONG),  Equivalent  slope  (EQSLP),  10/85  slope  (OSLP) , 
stream  length  (LEN),  years  of  daily  data  (DYRS),  and  Base  Flow 
Index  (BFI). 

It  is  evident  from  Table  2.1  that  for  extreme  value  analysis  of  low 
flows,  the  Gumbel  III  distribution  is  the  most  widely  accepted. 
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3.0  DATA  BASE 

3.1  General 

The  Southwestern-West  Central  regions  defined  by  the  Ministry  of 
Environment  were  used  as  the  study  area.  One  hundred  and  four 
hydrometric  stations  were  identified  in  the  region  for  this  analy- 
sis based  upon  the  period  of  record  and  other  station  characteris- 
tics. The  stations  considered  are  identified  in  Table  3.1.  Sixty- 
five  of  these  stations  were  used  while  ten  of  those  remaining  were 
set  aside  for  testing  of  results.  (Station  screening  and  selection 
is  discussed  in  Section  3.2.4.) 

Analyses  were  undertaken  to  produce  extreme  value  recurrence  low 
flows  and  tables  of  flow  exceedance  values  at  each  station.  The 
methodology  is  discussed  in  Section  3.2  and  results  are  summarized 
in  the  Low  Flow  Characteristics  in  Ontario*  report  (Appendix  D)  and 
on  Figure  3.1  of  this  report.  Station  physiographic  and  meteoro- 
logical characteristics  were  collected  and  are  discussed  in  Section 
3.3. 

3.2  Low  Flow  Characteristics 

3.2.1  Data  Base 

Existing  low  flow  characteristics  at  the  time  of  this  study  were 
several  years  out  of  date  and,  therefore,  in  conjunction  with  this 
study,  the  Low  Flow  Characteristic  map  was  updated  for  this 
region*. 


Cumming  Cockburn  Limited,  "Low  Flow  Characteristics  in  Ontario", 
Ministry  of  the  Environment,  1988 
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The  available  discharge  data  base  to  the  end  of  1985  was  obtained 
for  all  relevant  Water  Survey  of  Canada  monitoring  locations  (see 
discussion  in  Section  3,1).  Extreme  values  were  then  calculated 
for  various  recurrence  intervals  and  flow  durations  for  annual 
series  using  the  techniques  discussed  in  Section  3.2.5.  A  computer 
drawn  map  depicting  low  flow  characteristics  was  then  produced  in 
order  to  surmiarize  key  low  flow  statistics  for  Southwestern  West 
Central  regions  (see  Section  3.2.5  and  Figure  3.1). 

Data  analysis  and  screening  was  then  undertaken  in  order  to  assess 
the  usefulness  of  the  data  at  each  location. 

Stations  with  a  minimum  of  ten  years  of  record  and  stations  active 
within  the  last  five  years  were  considered  to  have  an  appropriate 
record  for  the  purpose  of  this  investigation.  These  criteria  were 
adopted  so  that  the  low  flow  characteristics  would  be  representa- 
tive of  recent  conditions. 

3.2.2  Extreme  Value  Analysis 

Average  extreme  low  flows  were  determined  and  extracted  for  the  1, 
3,  7,  15  and  30  day  durations  and  are  available  as  part  of  the 
background  files. 

An  extreme  value  analysis  was  undertaken  for  each  of  the  1,  3,  7, 
15  and  30-day  durations  for  each  of  the  stations. 

Previous  investigations  have  found  that  the  Gumbel  III  distribution 
has  resulted  in  the  best  fit  for  extreme  value  analysis  of  low 
flows  on  most  Canadian  rivers  (see  Section  2.0).  However,  for 
samples  with  large  negative  skewness,  the  3-Parameter  Log  Normal 
Distribution  has  proven  to  give  adequate  results  and  was-adopted 
for  these  stations  (Condie,  1983).  A  technical  discussion  of  the 
fitting  procedure  is  given  in  the  users  manual  and  briefly  in 
Appendix  A.l. 
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Parameter  estimation  for  fitting  of  distributions  proceeded  in 
order  of  maximum  likelihood,  smallest  observed  drought  and  mom- 
ents. The  procedure  utilized  was  identified  on  sunmary  tables 
together  with  the  results  of  the  analyses. 

The  results  of  the  extreme  value  analysis  for  each  station  are 
given  in  "The  Low  Flow  Characteristics  of  Ontario"  study  (Cuiming 
Cockburn  Limited,  1988)  for  various  recurrence  intervals  for  each 
duration.  The  format  is  illustrated  in  Table  3.2.  The  stations 
are  identified  by  the  Water  Survey  of  Canada  station  number.  The 
fitting  method  is  identified  by  a  3  letter  code,  MAX,  MOM,  SOD,  PLN 
which  stand  for;  method  of  maximum  likelihood,  method  of  moments, 
method  of  smallest  observed  drought,  and  the  three-parameter  log 
normal  distribution  respectively  (see  also  Appendix  A.l  for  addi- 
tional information  on  the  fitting  procedures). 

The  next  three  parameters  appearing  in  the  summary  table  are  gener- 
al statistics  of  the  data  sample  including  the  mean  n-day  duration 
flow  (m^/s),  Standard  deviation.  Skew  (G),  and  the  coefficient  of 
Variance  (C).  The  next  two  columns  of  the  tables  give  the  number 
of  years  available  (REC  (YRS))  and  the  minimum  average  low  flow  for 
the  particular  duration  (MIN). 

Finally,  the  estimated  flow  for  the  various  recurrence  intervals 
are  listed  for  each  station. 

Extreme  value  frequency  curves  were  plotted  for  each  location  and 
are  available  in  "The  Low  Flow  Characteristics  of  Ontario"  (Appen- 
dix D,  Cumming  Cockburn,  1988).  An  example  graph  is  given  in 
Figure  3.2  which  summarizes  the  curve  fitting  to  the  available  data 
for  the  1,  3,  7,  15  and  30  day  durations  for  the  Bronte  Creek  near 
Zimmerman,  Water  Survey  of  Canada  Station  Number  02HB011. 
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3.2.3  Flow  Duration  Analysis 

The  daily  discharge  data  was  extracted  for  each  station  and  analy- 
sed to  derive  flow  duration  curves.  The  flow  duration  curves  were 
summarized  both  numerically  and  graphically  for  each  station  ("The 
Low  Flow  Characteristics  of  Ontario"  (Appendix  D,  Gumming  Cockburn, 
1988)).  The  results  are  summarized  in  tabular  format  at  1%  inter- 
vals across  the  flow  duration  curve.  An  example  summary  table  is 
given  in  Table  3.3.  The  first  column  "PER"  refers  to  the  percent- 
age of  the  period  of  record  that  the  tabulated  flow  was  equalled  or 
exceeded  (all  flows  are  in  m^/s).  The  annual  flow  duration  inform- 
ation are  listed  in  column  2  for  the  percentage  summarized  in 
column  1.  Therefore,  the  largest  daily  flow  recorded  for  this  sta- 
tion up  to  1986  is  29.7  m^/s,  found  in  the  first  row  at  0  percent. 
More  significantly  for  this  study  are  the  low  flows  summarized  in 
the  latter  section  of  the  table  which  have  been  exceeded  90,  95  and 
100  percent  of  the  time  for  the  period  of  record. 

The  50  percent  value  (1.69  m^/s)  refers  to  the  median  daily  flow 
for  the  period  of  record  of  the  flow  series  and  can  be  compared  to 
the  mean  value  2.76  m^/s  summarized  in  the  last  row  of  the  table. 

3.2.4  Maps 

Selected  low  flow  characteristics  were  extracted  and  summarized  on 
a  map  for  this  region  (see  Figure  3.1).  The  station  locations  are 
identified  by  a  A  located  on  the  rivers  at  the  point  of  discharge 
measurement  and  the  selected  data  is  summarized  in  an  information 
box,  an  example  of  which  is  depicted  in  Figure  3.3. 

The  stations  are  identified  by  the  7  digit  Water  Survey  Station 
Number  and  this  is  followed  by  a  regulation  code.  The  "R"  indi- 
cates that  data  collected  at  the  station  is  affected  by  regulation, 
"N"  means  the  station  data  are  natural  or  non-regulated.  The  boxes 
on  the  left  from  the  top  refer  to  the  average  7-day  flow  (m^/s) 


EXAMPLE    TABLE    3-3 


SUMRY  TABLE  FROV  FLOi  DURATION  ANALYSIS  0H011  BROKTE  CREEK  )CMt  ZIMOUAN 

YEARS  OF  RECORD:        21  STATION  AREA:  235 

PER    ANNUAL       JANUARY    FEBRUARY     MARCH  APRIL  HAY          Jl*C          JULY         AUGUST  SEPTEiCER  XTOBER  NOVEMER  OECEICER 

0  29.700       22.700       23.000       28.900  29.700  21.500       10.300       12.900         9.080  12.900  15.700  10.500  17.500 

1  15.200         7.650       15.700       21.200  21.700  10.100         5.300         6.370         3.970  7.690  8.130  9.380  10.800 

2  12.600         6.460        14.300        17.600  19.100  8.670         4.470         4.050         3.430  5.780  6.740  8.340  8.210 


3 

10.700 

5.970 

13.200 

16.200 

17.100 

7.990 

3.990 

3.370 

3.060 

4.590 

6.160 

7.160 

7.700 

4 

9.710 

5.510 

12.400 

15.300 

15.500 

7.330 

3.670 

2.860 

2.520 

3.900 

5.590 

6.300 

7.190 

5 

8.880 

5.270 

11.300 

15.000 

14.900 

6.770 

3.420 

2.580 

2.410 

3.320 

5.130 

5.950 

6.800 

6 

8.240 

4.730 

9.630 

14.400 

14.000 

6.400 

3.340 

2.290 

2.290 

3.070 

4.670 

5.140 

6.440 

7 

7.650 

4.530 

8.410 

13.700 

13.400 

5.970 

3.230 

2.170 

2.140 

2.770 

4.500 

4.870 

6.170 

8 

7.080 

4.250 

7.890 

13.200 

12.700 

5.800 

3.110 

2.040 

2.000 

2.520 

4.020 

4.500 

6.000 

9 

6.560 

4.170 

7.080 

12.700 

12.100 

5.540 

3.060 

1.940 

1.940 

2.310 

3.820 

4.360 

5.790 

10 

6.210 

3.960 

6.570 

12.200 

11.600 

5.350 

2.990 

1.840 

1.760 

2.260 

3.510 

4.160 

5.460 

11 

5.860 

3.820 

5.950 

12.000 

11.200 

5.100 

2.900 

1.800 

1.700 

2.080 

3.310 

3.990 

5.320 

12 

5.520 

3.720 

5.660 

11.300 

10.700 

4.980 

2.810 

1.700 

1.590 

2.010 

3.050 

3.940 

5.210 

13 

5.230 

3.680 

5.380 

11.100 

10.500 

4.850 

2.750 

1.650 

1.540 

1.980 

2.330 

3.870 

5.1W 

14 

5.000 

3.540 

5.100 

10.900 

10.100 

4.690 

2.660 

1.610 

1.500 

1.910 

2.630 

3.790 

4.960 

15 

4.740 

3.400 

4.840 

10.500 

9.800 

4.620 

2.580 

1.580 

1.450 

1.830 

2.540 

3.740 

4.740 

16 

4.520 

3.310 

4.470 

10.500 

9.600 

4.470 

2.510 

1.530 

1.400 

1.750 

2.410 

3.710 

4.590 

17 

4.300 

3.250 

4.250 

10.300 

9.370 

4.370 

2.410 

1.490 

1.380 

1.710 

2.330 

3.650 

4.420 

18 

4.160 

3.170 

3.900 

10.1M 

9.150 

4.280 

2.380 

1.460 

1.360 

1.670 

2.220 

3.610 

4.320 

19 

4.000 

3.110 

3.680 

9.800 

8.960 

4.220 

2.250 

1.440 

1.310 

1.640 

2.160 

3.570 

4.250 

!0 

3.880 

3.030 

3.510 

9.670 

8.780 

4.120 

2.180 

1.380 

1.280 

1.600 

2.080 

3.540 

4.220 

!1 

3.740 

2.950 

3.450 

9.570 

8.690 

4.000 

2.130 

1.350 

1.240 

1.530 

2.020 

3.450 

4.190 

2 

3.650 

2.320 

3.370 

9.430 

8.500 

3.960 

2.1OT 

1.320 

1.220 

1.500 

1.940 

3.400 

4.110 

!3 

3.530 

2.890 

3.200 

9.300 

8.350 

3.910 

2.020 

1.310 

1.190 

1.470 

1.920 

3.340 

4.010 

!4 

3.430 

2.830 

3.140 

9.080 

8.160 

3.840 

1.980 

1.270 

1.180 

1.390 

1.880 

3.280 

3.940 

S 

3.340 

2.750 

2.990 

8.920 

8.040 

3.740 

1.950 

1.250 

1.160 

1.340 

1.820 

3.260 

3.890 

16 

3.2S0 

2.680 

2.S30 

8.800 

7.930 

3.660 

1.910 

1.230 

1.120 

1.310 

1.790 

3.200 

3.820 

7 

3.160 

2.610 

2.780 

8.600 

7.700 

3.620 

1.880 

1.200 

1.1W 

1.260 

1.770 

3.150 

3.770 

!8 

3.090 

2.550 

2.690 

8.500 

7.500 

3. SO 

1.830 

1.180 

1.080 

1.210 

1.700 

3.110 

3.740 

3.000  2.490  2.660  8.390  7.390  3.500  1.790  1.160  1.060  1.180  1.670  3.030  3.B50 

2.920  2.450  2.580  8.240  7.280  3.430  1.760  1.140  1.050  1.150  1.660  3.000  3.600 

2.830  2.410  2.550  8.070  7.160  3.390  1.710  1.110  1.030  1.100  1.640  2.950  3.590 

2.750  2.350  2.520  8.000  7.050  3.340  1.670  1.090  1.010  1.080  1.620  2.920  3.540 

2.670  2.270  2.460  7.800  6.930  3.280  1.660  1.080  1.000  1.060  1.610  2.890  3.460 

2.600  2.270  2.440  7.590  6.820  3.260  1.620  1.060  0.983  1.020  1.580  2.830  3.430 

2.520  2.270  2.380  7.370  6.680  3.230  1.590  1.040  0.964  0.994  1.550  2.760  3.400 

2.460  2.210  2.350  7.110  6.600  3.140  1.570  1.030  0.949  0.974  1.540  2.720  3.350 

2.410  2.190  2.290  7.000  6.460  3.090  1.530  1.010  0.940  0.951  1.500  2.660  3.280 

2.350  2.150  2.220  6.800  6.430  3.090  1.500  0.994  0.928  0.927  1.480  2.640  3.200 

2.290  2.120  2.120  6.650  6.400  3.020  1.480  0.974  0.914  0.912  1.450  2.600  3.150 


2.240 

2.100 

2.020 

6.430 

6.340 

2.970 

1.470 

0.966 

0.895 

0.898 

1.410 

2.560 

3.090 

2.170 

2.070 

1.930 

6.200 

6.210 

2.920 

1.440 

0.955 

0.882 

0.881 

1.370 

2.520 

3.060 

2.120 

2.040 

1.870 

6.040 

6.090 

2.890 

1.430 

0.934 

0.875 

0.871 

1.350 

2.470 

3.000 

2.060 

1.980 

1.700 

5.720 

6.030 
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0.867 

0.867 

1.310 

2.410 

2.920 
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1.980 

1.620 

5.660 

5.920 

2.810 

1.370 

0.911 

0.858 

0.852 

1.290 

2.380 

2.810 

1.940 

1.950 

1.620 

5.490 

5.850 

2.800 

1.350 

0.901 

0.841 

0.841 

1.270 

2.350 

2.750 

1.890 

1.900 

1.610 

5.380 

5.720 

2.760 

1.340 

0.886 

0.824 

0.827 

1.240 

2.300 

2.670 

1.840 

1.870 

1.590 

5.130 

5.610 

2.740 

1.320 

0.878 

0.816 

0.818 

1.210 

2.250 

2.650 

1.790 

1.840 

1.590 

5.060 

5.500 

2.700 

1.320 

0.850 

0.799 

0.812 

1.190 

2.210 

2.610 

1.75C 

1.810 

1.560 

4.990 

5.440 

2.660 

1.300 

0.835 

0.789 

0.8O4 

1.170 

2.150 

2.550 
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50 

1.690 

1.780 

1.550 

4.910 

5.340 

2.620 

1.290 

0.807 

0.779 

0.7S8 

1.150 

2.110 

2.520 

51 

1.650 

1.760 

1.530 

4.790 

5.250 

2.560 

1.280 

0.784 

0.783 

0.780 

1.130 

2.080 

2.500 

52 

1.620 

1.750 

1.490 

4.630 

5.210 

2.510 

1.250 

0.779 

0.755 

0.771 

1.100 

2.040 

2.460 

53 

1.590 

1.740 

1.460 

4.530 

5.130 

2.490 

1.250 

0.770 

0.744 

0.765 

1.080 

1.980 

2.430 

54 

1.550 

1.700 

1.430 

4.450 

5.030 

2.450 

1.230 

0.760 

0.738 

0.756 

1.060 

1.890 

2.390 

55 

1.520 

1.670 

1.420 

4.220 

4.330 

2.410 

1.230 

0.750 

0.728 

0.748 

1.030 

1.860 

2.370 

56 

1.490 

1.630 

1.420 

4.110 

4.S30 

2.390 

1.210 

0.733 

0.720 

0.739 

1.020 

1.830 

2.350 

57 

1.450 

1.610 

1.400 

3.990 

4.750 

2.350 

1.200 

0.728 

0.708 

0.731 

1.010 

1.760 

2.310 

58 

1.420 

1.590 

1.390 

3.940 

4.680 

2.340 

1.180 

0.722 

0.897 

0.725 

0.990 

1,730 

2.270 

59 

1.380 

1.560 

1.360 

3.780 

4.590 

2.310 

1.170 

0.711 

0.688 

0.716 

0.973 

1.680 

2.220 

60 

1.350 

1.540 

1.310 

3.710 

4.530 

2.270 

1.170 

0.708 

0.682 

0.711 

0.940 

1.S30 

2.170 

61 

1.310 

1.520 

1.270 

3.650 

4.470 

2.230 

1.150 

0.699 

0.677 

0.708 

0.920 

1.610 

2.150 

62 

1.290 

1.510 

1.250 

3.540 

4.400 

2.200 

1.140 

0.691 

0.668 

O.&Q 

0.900 

1.570 

2.120 

63 

1.250 

1.500 

1.250 

3.430 

4.300 

2.170 

1.120 

0.685 

0.657 

0.685 

0.895 

1.560 

2.  IX 

64 

1.220 

1.500 

1.200 

3.400 

4.250 

2.150 

1.110 

0.680 

0.651 

0.673 

0.883 

1.510 

2.040 

65 

1.190 

1.500 

1.190 

3.300 

4.220 

2.120 

1.100 

0.671 

0.640 

0.668 

0.869 

1.480 

1.980 

1.160  1.470  1.160  3. 260  4.160  2.080  1.090  0.663  0.630  0.S7  0.858  1.440  1.940 

1.130  1.420  1.130  3.260  4.080  2.080  1.080  0.657  0.623  0.646  0.850  1.400  1.890 

1.100  1.400  1.120  3.200  4.050  2.040  1.070  0.654  0.612  0.640  0.833  1.370  1.870 

1.070  1.360  1.100  3.140  3.990  2.020  1.050  0.648  0.606  0.634  0.827  1.350  1.850 


70 

1.060 

1.340 

1.080 

3.060 

3.940 

1.990 

1.060 

0.645 

0.600 

0.829 

0.809 

1.300 

1.810 

71 

1.020 

1.320 

1.060 

2.970 

3.850 

1.940 

1.CJ30 

0.631 

0.533 

0.623 

0.784 

1.260 

1.800 

72 

0.991 

1.310 

1.040 

2.890 

3.790 

1.930 

1.010 

0.623 

0.580 

0.623 

0.776 

1.240 

1.760 

73 

0.968 

1.300 

1.030 

2.800 

3.710 

1.900 

1.000 

0.617 

0.572 

0.612 

0.765 

1.190 

1.730 

74 

0.337 

1.280 

1.020 

2.760 

3.570 

1.880 

0.988 

0.612 

0.566 

0.602 

0.750 

1.180 

1.690 

75 

0.912 

1.220 

1.000 

2.730 

3.510 

1.840 

0.977 

0.602 

0.561 

0.598 

0.742 

1.150 

1.660 

76 

0.886 

1.190 

0.980 

2.660 

3.450 

1.810 

0.960 

0.592 

0.555 

0.589 

0.733 

1.120 

1.640 

77 

0.864 

1.150 

0.967 

2.550 

3.400 

1.800 

0.940 

0.583 

0.547 

0.581 

0.719 

1.080 

1.630 

78 

0.841 

1.120 

0.949 

2.450 

3.280 

1.770 

0.331 

0.568 

0.541 

0.573 

0.711 

1.060 

1.590 

79 

0.810 

1.070 

0.SJ4 

2.380 

3.230 

1.720 

0.923 

0.561 

0.538 

0.566 

0.705 

1.050 

1.560 

80 

0.782 

1.030 

0.920 

2.320 

3.180 

1.690 

0.905 

0.555 

0.532 

0.549 

0.692 

1.030 

1.540 

81 

0.760 

1.000 

0.860 

2.280 

3.140 

1.670 

0.888 

0.547 

0.521 

0.539 

0.679 

0.933 

1.500 

82 

0.736 

0.940 

0.830 

2.180 

3.090 

1.640 

0.867 

0.538 

0.515 

0.529 

0.668 

0.971 

1.460 

83 

0.716 

0.720 

0.800 

2.050 

3.030 

1.600 

0.861 

0.532 

0.500 

0.521 

0.66O 

0.960 

1.420 

84 

0.699 

0.600 

0.790 

1.950 

2.950 

1.570 

0.844 

0.523 

0.433 

0.513 

0.654 

0.915 

1.400 

85 

0.679 

0.572 

0.765 

1.850 

2.920 

1.550 

0.838 

0.513 

0.481 

0.499 

0.646 

0.887 

1.370 

86 

0.66O 

0.558 

0.700 

1.830 

2.860 

1.510 

0.810 

0.501 

0.476 

0.493 

0.637 

0.867 

1.330 

87 

0.643 

0.538 

0.570 

1.770 

2.780 

1.490 

0.784 

0.494 

0.470 

0.481 

0.633 

0.861 

1.300 

88 

0.626 

0.538 

0.634 

1.700 

2.740 

1.450 

0.776 

0.490 

0.464 

0.467 

0.623 

0.844 

1.300 

89 

0.612 

0.530 

0.595 

1.550 

2.670 

1.440 

0.759 

0.481 

0.459 

0.457 

0.614 

0.810 

1.300 

90 

0.595 

0.515 

0.540 

1.520 

2.610 

1.400 

0.736 

0.476 

0.450 

0.447 

0.603 

0.756 

1.260 

91 

0.570 

0.507 

0.510 

1.390 

2.520 

1.370 

0.708 

0.467 

0.442 

0.439 

0.595 

0.708 

1.210 

92 

0.549 

0.500 

0.498 

1.190 

2.510 

1.340 

0.695 

0.459 

0.433 

0.433 

0.567 

0.674 

1.180 

93 

0.533 

0.493 

0.490 

1.090 

2.430 

1.300 

0.681 

0.450 

0.428 

0.430 

0.561 

0.660 

1.150 

94 

0.513 

0.487 

0.480 

1.080 

2.360 

1.250 

0.657 

0.447 

0.413 

0.422 

0.544 

0.637 

1.020 

95 

0.496 

0.481 

0.453 

0.S34 

2.290 

1.200 

0.646 

0.439 

0.400 

0.418 

0.527 

0.620 

0.878 

96 

0.476 

0.479 

0.447 

0.900 

2.230 

1.170 

0.534 

0.433 

0.374 

0.399 

0.515 

0.600 

0.765 

97 

0.455 

0.464 

0.425 

0.738 

2.120 

1.120 

0.600 

0.408 

0.365 

0.382 

0.505 

0.583 

0.665 

98 

0.433 

0.456 

0.425 

0.640 

2.030 

1.060 

0.561 

0.394 

0.345 

0.374 

0.496 

0.564 

0.631 

99 

0.399 

0.311 

0.425 

0.59O 

1.780 

0.973 

0.515 

0.340 

0.328 

0.357 

0.470 

0.541 

0.614 

100 

0.286 

0.311 

0.419 

0.500 

1.610 

0.934 

0.433 

0.289 

0.286 

0.328 

0.433 

0.459 

0.600 

2.759         2.210  2.827         6.213         6.458         3.148  1.S33  1.124  1.010  1.243         1.706         2.469         3.092 
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with  a  recurrence  interval  of  2,  5,  10  and  20  years,  followed  by 
the  minimum  average  7-day  flow  and  the  period  of  record  for  the 
station.  The  values  shown  on  the  right  are  the  flows  (m^/s)  equal- 
led or  exceeded  for  the  available  period  of  record  5,  50,  75,  95, 
and  99  percent  of  the  time  followed  by  the  station  drainage  area  in 
km2. 

Station  names  are  listed  along  with  the  station  numbers  for  identi- 
fication purposes,  on  the  map  for  this  region. 

3.2.5  Data  Analysis  and  Screening 

The  available  data  for  each  station  were  reviewed  to  delete  sta- 
tions which  contained  30%  or  more  zero  low  flow  occurrences  for  the 
different  low  flow  durations.  The  stations  which  were  subsequently 
omitted  from  the  analysis  are  tabulated  and  discussed  in  Section 
3.2.6. 

A  draft  version  of  the  Low  Flow  Frequency  Analysis  program  was 
obtained  from  Environment  Canada  for  use  in  the  investigation 
(Pilon,  1987).  This  program  contains  a  set  of  data  analysis  and 
checking  modules  suitable  for  data  screening  and  analysis  (Pilon, 
1985).  However,  at  present  there  is  some  question  concerning  the 
methodology  for  checking  Independence  for  low  flow  data  series 
(personal  communication  Water  Resources  Branch,  Environment 
Canada)*. 

It  must  also  be  expected  that  a  certain  number  of  stations  which 
are  tested  will  randomly  "fail"  any  given  statistical  test.  There- 
fore, the  binomial  distribution  was  used  to  calculate  the  number  of 
expected  random  failures  over  the  entire  data  base  for  each  region 
in  order  to  assess  the  suitability  of  the  low  flow  data  on  a 
regional  basis.  A  tabulation  of  theoretical  and  actual  "failure" 
results  is  summarized  in  Table  3.4  and  discussed  in  Section  3.2.6. 


< 

{^ 

is 

cn 
o 

I- 

I- 
o 
< 
cr 

< 


i  3 


2  <r  < 

c^  o  - 

O    ^  Q 

w  o  S 

'*'  a 

9  2  £ 


O 
I 


C7>   a.    4   2    00 


ES 

ez' 

a 

g 

S.I 

UJ 

HI 

^ 

S 

O 

u" 

o 

o 

o 

9 

< 

1 

^  0.1  ^M 

3-6 

The  tests  were  also  re-evaluated  for  various  data  groupings  reflec- 
ting the  length  of  record  and  degree  of  regulation.  The  period  of 
record  was  subdivided  into  stations  with  less  than  20  years  of 
record  and  stations  with  greater  than  or  equal  to  20  years  of 
record.  The  effect  of  regulation  was  examined  by  further  subdivi- 
sion of  the  data  base  according  to  regulated  and  non-regulated  sta- 
tions (see  Table  3.4). 

3.2.6  Screening  Results 

Table  3.5  lists  the  stations  which  were  not  considered  to  possess  a 
suitable  low  flow  data  base  for  the  purposes  of  this  investiga- 
tion. Hence  these  stations  were  not  included  in  the  regionaliza- 
tion  analysis.  These  stations  were  excluded  for  the  following 
reasons: 

i)  Zero  Flows 

Prior  to  undertaking  the  extreme  value  analyses  all  data  sets 
containing  30%  or  more  zero  low  flows  for  all  durations  were 
identified  and  removed  from  the  data  base  (see  Table  3.3). 

ii)  Arithmetic  Overflow 

For  5  stations,  the  integral  of  the  probability  function 


F(X)  =  1  -  e.p  [-  l^f] 


*  The  trend  and  independence  test  assumed  a  lag  of  1  year  between  samples 
while  in  actual  fact  greater  than  1.5  year  or  less  than  0.5  year  lag 
between  samples  may  exist.  New  subroutines  were  created  to  exclude 
values  not  1  year  lag  from  testing  procedures,  re-analysing  the  data 
with  the  new  modules  is  beyond  the  scope  of  this  study. 
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where  e  is  the  lower  bound  parameter 
u  is  the  characteristic  drought 
a  is  the  shape  parameter 

could  not  be  obtained  since  the  function  causes  an  arithmetic  over- 
flow (i.e.  the  version  of  the  program  automatically  terminates  cal- 
culations). Hence,  extreme  value  analysis  could  not  be  undertaken 
using  such  data  without  program  modifications  which  were  beyond  the 
scope  of  this  investigation.  Therefore,  these  data  sets,  which  are 
identified  in  Table  3.5,  were  deleted  from  the  present  study. 
Additional  work  should  be  undertaken  to  modify  the  program  in  order 
to  obtain  suitable  curve  fitting  results. 

iii )  High  Outl iers 

The  analysis  of  stations  with  very  large  low  flows  (greater  than 
1000  m^/s)  could  not  be  completed  and  hence  were  identified  and 
removed  from  the  analysis  (e.g.  G2HA003  -  Niagara  River  at  Queens- 
ton).  Additional  research  investigations  are  required  to  identify 
reasons  for  program  termination  in  such  cases. 

iv)  Multiple  Stations 

It  was  determined  that  large  river  systems  may  have  several  stream- 
flow  gauges  at  several  locations  along  the  channel.  The  use  of  all 
such  highly  correlated  data  could  adversely  bias  the  development  of 
regional ization  methods.  Therefore,  only  representative  gauges  for 
such  streams  were  retained  for  further  analysis. 

v)  Heavily  Regulated 

Low  flows  on  heavily  regulated  streams  (i.e.  Thames  River  and  Grand 
River  systems)  are  affected  through  the  use  of  reservoirs  for  low 
flow  augmentation.  Since  the  main  focus  of  this  study  is  to  pro- 
duce a  method  for  investigating  low  flows  for  ungauged  watersheds. 


3-8 

stations  indicating  a  high  degree  of  regulation  were  excluded  from 
the  data  base.  Some  "regulated"  stations  were  retained  and  assign- 
ed a  regulatory  code  representing  degree  of  regulation.  Subsequent 
simple  correlation  analysis  confirmed  that  these  stations  could  be 
retained  in  the  data  set  due  to  insignificant  levels  of  correla- 
tion. 

vi)  Statistical  Tests 

Statistical  data  analysis  tests  were  undertaken  as  outlined  in  Sec- 
tion 3.2.5  and  described  in  more  detail  in  Appendix  A. 2.  The 
available  test  statistics  were  recently  made  available  as  part  of 
the  LFA  (Pilon  and  Jackson,  1987)  low  flow  package  and  their  appli- 
cation here  is  one  of  the  first  in  an  investigation  of  low  flow 
characteristics. 

In  general,  it  was  found  that  a  significant  number  of  stations 
"failed"  the  non-parametric  tests.  Therefore,  taken  over  the 
entire  data  base,  application  of  these  tests  has  indicated  that  the 
available  data  base  of  extreme  low  flows  may  exhibit  some  trend  and 
dependence  with  some  possibility  of  non-random  characteristics. 

The  data  were  further  analysed  by  subdivision  of  the  available  data 
set  according  to  length  of  record  (i.e.  >20  years  and  <20  years) 
and  according  to  regulation  code.  However,  it  was  found  that 
neither  the  length  of  record,  nor  the  possible  effects  of  regula- 
tion could  account  for  the  conclusions  of  the  test  results.  One 
explanation  could  be  that  the  available  record  lengths  are  too 
short  to  permit  reasonable  application  and  interpretation  of  these 
non-parametric  test  results.  A  stronger  possibility  is  that  the 
available  low  flow  data  sets  do  exhibit  trend  and  non  random  char- 
acteristics, which  could  possibly  be  attributed  to  slow  cyclic 
change  in  groundwater  levels  or  to  climatic  trends.  Additional 
testing  was  beyond  the  scope  of  the  current  investigations. 
However,  further  studies  are  recommended  since  these  results  may 
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call  into  question  the  basic  assumptions  underlying  application  of 
the  extreme  value  analysis  technique  for  analysis  of  low  flow 
characteristics. 

Stations  not  passing  the  statistical  screening  were  not  removed 
from  the  analysis  as  there  would  be  too  few  stations  remaining. 

3.2.7  Summary  of  Data  Analysis  and  Low  Flow  Characteristics 

Application  of  various  non-parametric  tests  was  undertaken  for  the 
available  data  base  for  the  first  time  in  Ontario.*  The  test 
results  indicate  that  the  available  data  base  of  extreme  low  flows 
may  exhibit  some  trend  and  dependence  with  some  possibility  of 
non-random  characteristics.  Previous  wide-spread  application  of 
the  tests  utilized  have  not  been  found  in  the  literature  for  low 
flows.  The  average  length  of  record  for  the  stations  analysed  in 
this  region  is  26  years.  It  is  possible  that  the  available  record 
lengths  are  too  short  to  permit  reasonable  application  and  inter- 
pretation of  these  non-parametric  test  results.  Another  possibil- 
ity, which  should  be  investigated  in  more  detail,  is  that  the 
available  low  flow  data  sets  do  exhibit  trend  and  non-random  char- 
acteristics. The  latter  could  possibly  be  attributed  to  slow  cyc- 
lic changes  in  groundwater  levels  or  climatic  trends.  This  should 
be  investigated  in  more  detail  since  the  underlying  assumptions  of 
the  extreme  value  and  subsequent  regression  analyses  are  called 
into  question. 

The  Gumbel  Extreme  Value  Distribution  was  found  to  adequately  fit 
the  available  low  flow  series  for  various  low  flow  durations.  Data 
for  some  stations  could  not  be  fit  due  to  either  a  large  number  of 
zero  flows,  numerical  computational  problems,  or  very  large  low 


*  The  analysis  was  done  in  conjunction  with  "Low  Flow  Characteristics  in 
Ontario  Study",  Cumming  Cockburn  Limited,  1988 
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flows  recorded  at  some  locations.  Additional  research  investiga- 
tions are  required  to  develop  suitable  low  flow  analysis  techniques 
for  such  locations. 

Flow  duration  analysis  were  undertaken  on  a  total  of  just  over  100 
stations,  both  on  an  annual  and  monthly  basis.  Extreme  value  anal- 
ysis were  also  undertaken  on  an  annual  basis  for  just  over  100  sta- 
tions. A  total  of  65  stations  were  retained  for  analysis  and  10 
for  additional  testing  of  results. 

The  data  analyses  were  undertaken  for  both  regulated  and  unregula- 
ted data  series.  Therefore  care  should  be  taken  in  comparison  and 
interpretation  of  results,  notably  for  data  series  which  may 
include  the  effects  of  regulation. 

It  was  found  that  low  flows  are  consistently  lower  in  the  South- 
western/West Central  regions  than  the  Provincial  average.  Mean  low 
flow  averages  for  the  Province  range  from  10  to  15  m^/s  while  they 
were  found  to  range  from  2  to  3  m^/s  for  the  Southwestern  West 
Central  Region  for  1  to  30  day  averages.  The  mean  7-day  low  flow 
for  all  stations  analysed  for  this  region  was  found  to  be  2.44  m^/s 
and  the  mean  unit  area  7Q20  low  flow  is  0.70  l/s/km^  with  a  stand- 
ard deviation  of  0.86.  The  regional  differences  of  mean  7  day  low 
flow  values  may  be  attributed  to  a  function  of  the  hydrometric  net- 
work regional  differences  in  watershed  areas. 

Figure  3.1  summarizes  the  following  low  flow  characteristics: 

-  7  day  extreme  values  for  the  2,  5,  10  and  20  year  recurrence 
intervals 

-  flows  which  were  equalled  or  exceeded  over  the  available 
period  5,  50,  75,  95  and  99  percent  of  the  time 

Data  analysis  and  management  technique  are  now  available  which 
would  allow  efficiently  updating  the  present  analyses  on  a  frequent 
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basis.  In  our  opinion,  the  low  flow  analyses  should  be  updated 
every  five  years  in  order  to  provide  reasonably  accurate 
information  for  investigations  requiring  low  flow  information. 

3.3    Physiographic  and  Hydrometeoroloqic  Data 

3.3.1  Criteria 

Several  criteria  were  used  to  identify  physiographic  and  climatic 
parameters  which  might  be  suitable  for  use  in  regionalizing  low 
flow  characteristics.  These  are  discussed  as  follows: 

i)  Statistical  Significance 

When  undertaking  a  multivariate  analysis,  the  variables  chosen  must 
make  a  contribution  to  explaining  the  variance  of  the  low  flows. 
The  significant  experience  in  undertaking  similar  investigations 
(discussed  in  Section  2.0)  was  used  to  identify  parameters  which 
have  proven  to  be  statistically  significant  in  predicting  low 
flows. 

ii)  Physical  Reality 

Wherever  possible,  variables  should  be  selected  based  on  hydrologic 
significance.  That  is,  the  parameters  should  have  some  physical 
meaning  with  regard  to  estimate  of  low  flows. 

iii)  Reliability  of  Computation 

It  is  preferable  to  select  parameters  which  can  easily  be  computed 
in  a  reliable  manner  by  users  who  may  not  be  familiar  with  regress- 
ion procedures  or  the  details  of  the  statistical  concepts.  There- 
fore from  a  practical  point  of  view,  it  was  desirable  to  make  the 
parameter  estimation  procedure  as  uncomplicated  as  possible  in 
order  to  minimize  computational  errors  when  applying  the  technique. 
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3.3.2  Parameters 

The  parameters  selected  for  use  in  this  study  are  listed  as 
fol lows: 

Hydrometeoroloqic  Data  Symbol 
Index  of  mean  annual  precipitation  at  gauge 

location  (mm)  MAP 

Index  of  mean  annual  snowfall  at  gauge  location  (cm)  MAS 

Index  of  mean  annual  runoff  at  gauge  location  (mm)  MAR 

Physiographic  Data  Symbol 

Drainage  area  (km^)  DA 

Index  of  area  controlled  by  lakes  and  swamps  ACLS 

Latitude  of  gauge  location  LAT 

Longitude  of  gauge  location  LONG 

Length  of  main  channel  (km)  LEN 

10/85  channel  slope  (m/m)  SLR 

Base  Flow  Index  (dimensionless  ratio)  BFI 
Regulation  Index  (0  -  non-regulated,  3  -  slightly 
regulated,  6  -  moderately  regulated,  9  -  heavily 

regulated)  REG 
Regional  Code  (key  used  in  identifying  homogeneous 
regions  1  -  Lake  Huron  Region 

2  -  West  Lake  Erie 

3  -  East  Lake  Erie  RG 

Detailed  parameter  definition  and  methodology  for  derivation  is 
discussed  in  the  following  sections. 

DA  (km^) 

The  watershed  drainage  area  was  obtained  from  records  published  by 
the  Water  Survey  of  Canada.  However,  in  a  few  cases,  it  was 
necessary  to  modify  published  drainage  areas  for  a  number  of 
reasons  (e.g.  diversions,  more  accurate  mapping,  etc.). 
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ACLS 

An  index  representing  the  percentage  of  the  drainage  area  controll- 
ed by  lakes  and  swamps  (ACLS)  was  determined  using  available  topo- 
graphic mapping  for  all  basins  (Scale  1:50,000).  A  lake  or  swamp 
had  to  have  a  surface  area  equal  to  at  least  1%  of  the  drainage 
area  to  the  outlet  of  the  lake  or  swamp  in  order  to  be  considered 
to  have  a  significant  effect  on  low  flows. 


LAT 


The  station  latitude  expressed  in  degrees  and  decimal  format  was 
obtained  from  Water  Survey  of  Canada  publications. 


LONG 


The  longitude  of  the  hydrometric  station  expressed  in  degrees  and 
decimal  format  was  obtained  from  records  published  by  the  Water 
Survey  of  Canada. 

SLOPE  (m/m) 


Channel  slopes  were  determined  by  creating  an  elevation  distance 
table  which  was  measured  along  the  main  stream  of  each  river  system 
using  1:50,000  scale  mapping. 

The  10/85  slope  (SLP)  recommended  by  the  U.S.  Geological  Survey  was 
calculated  by  taking  the  elevation  difference  at  the  10%  and  85% 
points  along  the  channel  length. 

LEN  (km) 

The  length  of  the  main  channel  was  determined  from  1:50,000  scale 
NTS  mapping. 
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Base  Flow  Index  (BFI)  (Dimensionless) 

This  parameter  is  an  indicator  of  the  hydrogeological  effects  of 
the  drainage  basin  soil  and  geology  type  and  also  the  retention 
characteristics  (primarily  due  to  lakes  and  swamps)  of  the  drainage 
basin.  BFI  is  defined  as: 

BFI  =  Total  volume  of  base  flow/total  volume  of  runoff 

The  median  values  calculated  for  all  Ontario  gauging  stations  hav- 
ing at  least  2  years  of  continuous  daily  discharge  data  were  plott- 
ed at  the  corresponding  drainage  basin  centroids  using  1:600,000 
scale  base  maps  for  Southern  Ontario  and  isolines  drawn.  The  cent- 
roids were  located  by  eye  after  delineating  the  drainage  area.  The 
isolines  map  was  prepared  to  help  provide  estimates  of  BFI  for  un- 
gauged  basins  when  applying  the  regression  equations  (Regional 
Flood  Frequency,  Moin  and  Shaw,  1986).  All  estimates  of  BFI  from 
the  isoline  maps  must  be  made  by  first  locating  the  basin  centroid 
and  then  projecting  this  point  to  the  closest  point  on  the  main 
channel.  The  BFI  is  then  interpolated  from  the  isolines  at  this 
location  on  the  main  channel.  A  better  BFI  estimate  will  be  ob- 
tained for  large  basins  and  in  areas  where  the  isolines  are  very 
close  together,  if  an  average  value  of  BFI,  weighted  by  the  area 
between  isolines,  is  taken  over  the  entire  drainage  basin.  Figure 
3.4  reproduces  the  BFI  distribution  for  the  Southwestern  West  Cen- 
tral region. 


RG 


Identifying  regions  by  indexing  to  a  common  number  for  sub-region 
analysis.  Cluster  analysis  of  common  regression  residuals 
indicated  3  possible  homogeneous  regions  identified  as: 

1.  Lake  Huron  Region 

2.  West  Lake  Erie 

3.  East  Lake  Erie 
(see  Figure  4.1) 
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REG  (code) 

Regulation  code  (i.e.  an  assessment  of  the  degree  of  regulation  a 
stream  is  determined  to  have): 

0  -  non-regulated 

3  -  some  regulation 

6  -  moderate  regulation 

9  -  heavy  regulation 

NOTE:   Heavily  regulated  stations  were  not  included  in  the  analy- 
sis and  simple  correlation  analysis  indicated  insignificant 
levels  of  correlation  with  the  remaining  regulation  codes. 

MAP  (mm) 

An  index  of  mean  annual  precipitation  was  developed  with  reference 
to  available  publications  (MNR,  1984,  Environment  Canada,  1978). 
The  index  of  mean  annual  precipitation  is  interpolated  for  each 
watershed  used  in  the  analysis  at  the  gauge  location. 

MAS  (cm) 

An  index  of  mean  annual  snowfall  was  obtained  from  available 
information  published  by  the  Ministry  of  Natural  Resources  (MNR, 
1984).  An  index  map  was  used  to  derive  the  snowfall  index  for  each 
hydrometric  gauge.  This  index  represents  total  annual  snowfall  for 
each  watershed.  The  index  was  determined  at  the  gauge  location. 

MAR  (mm) 

The  index  of  mean  annual  runoff  is  expressed  as  a  depth  of  water 
averaged  over  the  drainage  basin  area.  The  isolines  of  runoff  were 
obtained  from  existing  information  published  by  Environment  Canada 
and  adopted  by  the  Ministry  of  Natural  Resources  (Sangal  and 
Kallio,  1977,  MNR,  1984).  The  mean  annual  runoff  index  was  derived 
from  for  each  hydrometric  station  at  the  gauge  location. 
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3.3.3  Summary  of  Physiographic  and  Hydrometeorologic  Data 

Table  3.6  summarizes  the  parameters  which  were  determined  for  each 
station  remaining  after  screening  (see  Section  3.2.5).  Table  3.7 
summarizes  the  mean,  range  and  other  simple  statistics  of  the  data 
base.  The  selected  data  contains  stations  with  drainage  area  vary- 
ing from  14  to  4000  km^,  with  BFI  which  range  from  .1  to  .8, 
lengths  ranging  from  6  to  190  km  and  slopes  ranging  from  .0003  to 
.0075  m/m. 

A  matrix  of  simple  correlation  coefficients  is  presented  in  Table 
3.8.  Significant  correlation  coefficient  has  a  magnitude  equal  or 
larger  than  0.24.  It  was  found  that  a  few  parameters  exhibit 
significant  intercorrelation  (e.g.  REG/LONG;  MAP/LAT;  MAP/MAS; 
DA/LEN;  DA/YEARS). 

Also,  latitude  was  found  to  be  highly  correlated  with  the  climato- 
logical  parameters  and  longitude  was  negatively  correlated  with  the 
cl imatological  parameters.  This  is  a  function  of  the  hydrometric 
network  and  does  not  represent  a  significant  hydrologic  finding  in 
regard  to  estimation  of  low  flows. 

As  one  would  intuitively  expect,  drainage  area  was  found  to  be 
highly  correlated  with  length.  This  finding  is  consistent  with 
other  investigations.  Generally  speaking,  equations  with  drainage 
area  should  not  likely  also  contain  length  as  an  important  para- 
meter, unless  the  regression  coefficients  are  found  to  be  very 
stable. 

While  not  shown  in  the  simple  correlation  matrix,  it  was  also  found 
that  the  low  flows  were  most  highly  correlated  with  drainage  area 
(0.85)  followed  by  length  and  BFI.  Low  flows  were  also  found  to  be 
negatively  correlated  with  slope,  a  fact  which  may  be  attributed  to 
large  base  flows  generated  from  large  watersheds  with  smaller 
slopes. 


Table  3.6 
■  SUMMARY    OF  INDEPENDENT  PARAMETERS 

R6    HAP    HAS    HAR    DA     ACLS   LAT     L0N6    LEN    SLP    YEARS  OF   BFI     REG 

RECORD 


(code) 

(M) 

(CI) 

(■■) 

(ka'2) 

(t) 

(decree) 

(degree) 

(ki)   (■/■) 

(years) 

(index) 

(code 

02FA001 

1 

910 

300 

463 

927.00 

98 

44.669 

81.253 

82.40  .0008900 

28 

.620 

OeFAOOS 

1 

890 

350 

325 

50.50 

80 

45.036 

81.336 

8.00  .0035700 

10 

.500 

02FBO07 

2 

1000 

300 

503 

181.00 

58 

44.700 

80.950 

27.60  .0024500 

47 

.560 

OEFBOIO 

2 

9W 

300 

481 

293.00 

49 

44.569 

80.650 

33.30  .0034800 

28 

.550 

02FC001 

1 

900 

300 

449 

3960.00 

1 

44.456 

81.327 

190.50  ,0017300 

71 

.640 

02FC011 

2 

1000 

280 

404 

163.00 

1 

44.113 

81.020 

29.50  .0020900 

32 

.510 

02FC0I2 

2 

1000 

280 

473 

635.00 

5 

44.099 

80.986 

97.50  .0023100 

14 

.420 

02FC013 

2 

920 

110 

255 

262.00 

30 

44.318 

81.246 

40.00  .0013200 

14 

.650 

02FC015 

2 

900 

290 

516 

663.00 

91 

44.269 

81.269 

83.00  .0014400 

14 

.620 

02FC016 

2 

1020 

270 

375 

329.00 

36 

44.186 

80.787 

52.30  .0027000 

9 

.660 

02FE002 

1 

960 

260 

421 

1630.00 

5 

43.887 

81.327 

76.90  .0012900 

32 

.410 

02FE003 

2 

920 

250 

386 

77.70 

9 

43.727 

80.972 

20.20  .0019800 

33 

.280 

02FE005 

1 

960 

260 

425 

528.00 

12 

43.916 

81.264 

32.40  .0013200 

32 

.460 

02FE007 

1 

950 

260 

464 

326.00 

1 

43.854 

81.250 

50.40  .0017100 

IB 

.400 

02FE008 

1 

9*0 

260 

432 

648.00 

2 

43.813 

81.307 

81.70  .0010000 

18 

.330 

02FE009 

1 

900 

250 

485 

373.00 

1 

43.684 

81.544 

46.50  .0014800 

18 

.320 

02FF002 

1 

890 

190 

358 

865.00 

1 

43.072 

81.660 

100.90  .0010600 

38 

.370 

02FF007 

1 

900 

240 

397 

466.00 

1 

43.300 

81.800 

52.50  .0019200 

19 

.330 

026A010 

3 

900 

170 

334 

1030.00 

18 

43.191 

80.455 

145.50  .0009600 

42 

.490 

026A018 

2 

890 

200 

342 

552.00 

1 

43.378 

80.711 

68.60  .0010200 

35 

.300 

026A024 

3 

840 

205 

350 

59.60 

65 

43.469 

80.519 

15.00  .0028700 

25 

.440 

026A02? 

3 

820 

180 

336 

236.00 

20 

43.548 

80.183 

42.00  .0043100 

23 

.560 

02BA030 

3 

880 

190 

137 

49.70 

1 

43.550 

80.250 

15.10  .0052800 

19 

.360 

02eA031 

3 

820 

180 

419 

44.50 

1 

43.576 

80.110 

6.10  .0043900 

20 

.730 

026A033 

3 

840 

200 

344 

64.80 

10 

43.663 

80.255 

11.40  .0034900 

29 

.460 

026A035 

2 

870 

280 

370 

27.70 

1 

43.651 

80.572 

16.00  .0074700 

14 

.430 

026A036 

2 

890 

240 

418 

17.90 

1 

43.670 

80.597 

6.50  .0059700 

14 

.200 

026A037 

3 

860 

200 

445 

25.10 

1 

43.439 

80.476 

6.30  .0049100 

13 

.320 

026A038 

2 

900 

200 

399 

326.00 

1 

43.484 

80.835 

27.00  .0016600 

13 

.190 

02eA039 

2 

840 

200 

417 

272.00 

1 

43.783 

80.639 

37.00  .0025000 

13 

.200 

026AO40 

3 

830 

200 

398 

167.00 

4 

43.639 

80.270 

22.50  .0037600 

12 

.390 

926B006 . 

3 

900 

150 

379 

150.00 

! 

43.175 

80.553 

42.50  .0016300 

32 

.480 

02SBOO7 

3 

900 

140 

309 

360.00 

20 

43.147 

80.155 

47.90  .0022100 

21 

.420 

D26B00e 

3 

900 

150 

345 

383.00 

1 

43.126 

80.385 

26.40  .0013000 

24 

.550 

D2eB010 

3 

900 

170 

301 

171.00 

1 

43.034 

79.950 

60.00  .0008300 

24 

.360 

m> 


HAS 


Table  3.6  (continued) 
SUMMARY    OF   INOEPENC)ENT    PARAMETERS 
HAR         DA  ACL5        LAT  L0N6  LEN 


SLP         YEARS  OF        BFI  RES 

RECORD 


(code) 

(M) 

(CI) 

(N) 

(kr2) 

(X) 

(degree) 

(degree) 

(ki) 

(■/I) 

(years) 

(index) 

(code) 

026C002      2 

890 

120 

369 

329.00 

42.778 

81.214 

45.50 

.0016800 

18 

.260 

026C006      3 

900 

140 

328 

363.00 

42.838 

80.510 

51.60 

.0011900 

30 

.690 

026C007      3 

890 

120 

335 

591.00 

42.685 

80.538 

85.00 

.0010600 

30 

.750 

028C008      3 

m 

130 

374 

134.00 

42.823 

80.289 

13.00 

.0025500 

28 

.730 

026C010      2 

900 

140 

354 

342.00 

42.861 

80.719 

40.50 

.0012000 

25 

.520 

026C0i2      3 

890 

140 

419 

51.30 

42.854 

80.328 

12.00 

.0017200 

21 

.800 

026C013      3 

900 

110 

401 

75.90 

42.628 

80.472 

19.00 

.0028000 

20 

.480 

026C015      2 

900 

120 

438 

105.00 

42.756 

80.807 

25.50 

.0017500 

22 

.670 

026C017      3 

900 

150 

390 

93.20 

42.965 

80.543 

10.50 

.0014300 

21 

.500 

026C018      2 

890 

110 

395 

287.00 

42.746 

81.044 

31.50 

.0013900 

21 

.280 

026C021       3 

900 

110 

250 

68.40 

42.674 

80.600 

22.50 

.0023600 

19 

.750 

026C026      3 

900 

100 

248 

676.00 

42.708 

80.839 

82.00 

.0007400 

10 

.620 

026D001       2 

900 

160 

336 

1340.00 

10 

42.973 

81.210 

88.80 

.0012500 

70 

.510 

026D004      2 

900 

200 

369 

306.00 

43.059 

80.988 

37.90 

.0021300 

37 

.400 

026D009      2 

900 

222 

350 

140.00 

100 

43.264 

81.080 

32.00 

,0017000 

34 

.440 

026D011       2 

890 

130 

315 

93.20 

43.128 

80.768 

12.80 

.0018300 

32 

.360 

026D012      2 

890 

180 

347 

254.00 

100 

43.144 

80.751 

36.80 

.0015400 

33 

.410 

026D015      2 

900 

210 

386 

1340.00 

10 

43.149 

81.192 

66.50 

.0010700 

32 

.370 

026D018      2 

890 

230 

437 

144.00 

43.345 

81.116 

24.60 

.0012300 

21 

.340 

026D020      2 

900 

180 

315 

108.00 

42.994 

81.117 

32.50 

.0025200 

20 

.360 

026E005      2 

890 

150 

316 

146.00 

42.934 

81.351 

37.20 

.0009500 

20 

.290 

0266002      2 

840 

140 

320 

730.00 

42.830 

81.853 

79.00 

.0008400 

37 

.420 

0266004       1 

780 

140 

288 

609.00 

42.764 

82.342 

89.50 

.0003600 

19 

.270 

0266005      1 

890 

160 

356 

172.00 

42.960 

81.628 

81.00 

.0008100 

19 

.470 

0266006      1 

800 

140 

299 

267.00 

42.906 

82.120 

36.50 

.0009400 

19 

.200 

0266007      2 

800 

120 

312 

1240.00 

42.594 

82.109 

57.50 

.0003400 

16 

.370 

026H001       2 

810 

90 

413 

14.20 

11 

42.046 

83.575 

8.10 

.0044500 

11 

.390 

026H002      2 

810 

110 

286 

125.00 

42.232 

82.094 

20.00 

.0009900 

14 

.150 

026H003      2 

790 

105 

200 

159.00 

25 

42.159 

83.018 

22.00 

.0004500 

9 

.100 

02HA006       3 

800 

120 

302 

293.00 

1 

43.134 

79.383 

80.00 

.0005500 

28 

.150 

Table  3.7 
Siiple  Statistics  of  Independent  Paraieters 


Paraieter 

RG  (code) 

HAP  (u) 

HAS  (CI) 

HAR  (■■) 

DA  (kt'^a) 

ACLS  (percent) 

LAT  (degrees) 

L0N6  (degrees) 

LEN  (ki) 

SLP  10/85  (■/■) 

YEARS  OF  RECORD  (years) 

6FI  (index) 

REB  (code) 


Mean   Std  Dev  Skewness  S.E.  SkeR    Range  Hiniiua  Haxiiui 


2. IS  .72  -.19 

888.62  50.37  .13 

189.72  64.09  Ab 

367.74  71.28  -.« 

414.00  571.57  4.15 

14.05  26.86  2.32 

43.34  .63  .54 

80.99  .70  1.11 

45.89  34.55  1.63 

.00  .00  1.69 

24.37  12.00  1.82 

.44  .16  .22 

2.75  2.16  .69 


.30 

2.00 

1 

3 

65 

,30 

240.00 

780 

1020 

65 

,30 

260.00 

90 

350 

65 

,30 

379.00 

137 

516 

65 

.30 

3945.80 

14.20 

3960.00 

65 

,30 

99.00 

1 

100 

65 

,30 

2.99 

42.046 

45.036 

65 

,30 

4.19 

79.383 

83.575 

65 

,30 

184.40 

6.10 

190,50 

65 

,30  7 

.130E-03 

.0003400 

.0074700 

65 

,30 

62.00 

9 

71 

65 

30 

.70 

.100 

.800 

65 

,30 

5.00 

1 

6 

65 

Table  3.8 
Correlation  Hatrix  of  Independent  Paraieters 


RG 

HAP 

HAS 

HAR 

DA 

ACLS 

LAT 

LONG 

LEN 

SLP 

RECORD 
LENGTH 

BFI 

RES 

(code) 

(■■) 

(CI) 

(■■) 

(kr2) 

(X) 

(degree) 

(degree) 

(ki) 

(■/I) 

(years) 

(index) 

(code) 

RG 

1 

-0.2013 

-0.4758 

-0.3159 

-0.3686 

-0.1312 

-0.3204 

-0.645 

-0.3813 

0,2605 

-0.0989 

0.312 

0.3253 

HflP 

-0.2013 

1 

0.5494 

0.4402 

0.1076 

0.1472 

0.5414 

-0.1982 

0.1465 

-0.0305 

0.1935 

0.3237 

0.142 

HAS 

-0.4758 

0.5494 

1 

0.5727 

0.2451 

0.4229 

0.8763 

-0.0534 

0.1792 

0.2707 

0.1821 

0.0672 

-0.0721 

NAR 

-0.3159 

0.4402 

0.5727 

1 

0.1896 

0.181 

0.4763 

-0.0157 

0.1377 

0.0586 

0.1592 

0.2004 

0.0009 

DA 

-0.3686 

0.1076 

0.2451 

0.1896 

1 

-0.0509 

0.2131 

0.1872 

0.8247 

-0.2831 

0.5843 

0.1057 

-0.137 

ACLS 

-0.1312 

0.1472 

0.4229 

0.181 

-0.0509 

1 

0.4596 

0.0308 

-0.0468 

0.0569 

0.0364 

0.1813 

-0.1453 

LAT 

-0.3204 

0.5414 

0.8763 

0.4763 

0.2131 

0.4596 

1 

-0.2246 

0.1635 

0.2661 

0.1408 

0.2208 

-0.0731 

LONG 

-0.645 

-0.1982 

-0.0534 

-0.0157 

0.1872 

0.0308 

-0.2246 

1 

0.1475 

-0.2637 

-0.1163 

-0.352 

-0.3443 

LEN 

-0.3813 

0.1465 

0.1792 

0.1377 

0.8247 

-0.0468 

0.1635 

0.1475 

1 

-0.4416 

0.4544 

0.0512 

-0.1748 

SLP 

0.2605 

-0.0305 

0.2707 

0.0586 

-0.2831 

0.0569 

0.2661 

-0.2637 

-0.4416 

1 

-0.1779 

0.1869 

0.056 

YEARS 

-0.0989 

0.1935 

0.1821 

0.1592 

0.5843 

0.0364 

0.1408 

-0.1163 

0.4544 

-0.1779 

1 

0.1857 

0.0375 

BFI 

0.312 

0.3237 

0.0672 

0.2004 

0.1057 

0.1813 

0.2208 

-0.352 

0.0512 

0.1869 

0.1857 

1 

0.2016 

RE6 

0.3253 

0.142 

-0.0721 

0.0009 

-0.137 

-0.1453 

-0.0731 

-0.3443 

-0.1748 

0.056 

0.0375 

0.2016 

1 
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4.0  REGIONALIZATION  OF  LOW  FLOW  CHARACTERISTICS 

4.1  General 

Four  alternative  methods  of  regionalizing  extreme  value  low  flows 
and  flow  duration  values  at  the  95%  level  (O95)  were  examined. 
They  are;  Multiple  Linear  Regression  (Section  4.2),  Index  method 
(Section  4.3),  Area  Proration  mapping  of  unit  flow  values  (Section 
4.4),  and  Isolines  of  unit  low  flows  (Section  4.5).  Emphasis  was 
placed  on  determining  7O20  and  O95  as  these  flows  were  identified 
to  be  the  key  low  flow  statistic  required  byr  M.O.E. 

4.2  Regression 


4.2.1  General 


Multiple  linear  regression  equations  take  the  following  general 
form: 

Y  =  ao  +  aiYi  +  a2Y2  +  +  apYn  (4.1) 

where        Y  =  the  dependent  variable  (e.g.  702o) 

Yi,Y2...Yn  =  the  independent  variables  (e.g.  physiographic 
and  hydrometeorologic  watereshed  character- 
istics) 
ao,ai...an  =  regression  coefficients 

In  order  to  obtain  a  more  accurate  regression  equation,  it  is  some- 
times necessary  to  transform  the  data  by  taking  logarithms,  square 
roots,  cubes,  etc.  The  transformations  considered  for  this  invest- 
igation are  discussed  in  Section  4.2.3. 

A  number  of  regression  procedures  are  available,  including  develop- 
ment of  all  possible  equations,  forward  selection,  backward  selec- 
tion, stagewise  regression  and  stepwise  regression.  The  stepwise 
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regression  procedure  is  generally  recormiended  for  use  in  practical 
applications  (Draper  and  Smith,  1981)  and  was,  therefore,  adopted 
for  the  purposes  of  this  investigation. 


4.2.2  Methodology 


The  regression  equations  were  developed  utilizing  the  stepwise  mul- 
tiple linear  regression  procedure  available  in  the  Statistical 
Package  for  the  Social  Sciences  (SPSS)  (Nie,  1975).  More  specific- 
ally, the  regression  subprogram  has  been  used  and  has  the  following 
special  features: 

1.  Various  procedures  for  selecting  variables,  including  forward 
selection,  backward  elimination  and  the  stepwise  selection 
feature,  which  is  really  a  combination  of  backward  and  for- 
ward procedures,  was  adopted  for  this  investigation  as  it  is 
the  most  comnonly  used  procedure. 

2.  Variable  selection  -  all  independent  variables  can  be  stored 
and  then  only  those  variables  desired  for  a  particular  analy- 
sis called  up  and  used  according  to  the  desired  form  of  the 
equation 

3.  Combination  of  variables  -  variable  transformation  and  new 
variables  may  be  computed  from  existing  variables 

4.  Transformations  -  the  variables  may  be  transformed  (e.g. 
square  root,  logarithmic,  reciprocal,  etc.)  in  order  to  more 
nearly  linearize  the  relationships 

5.  Calculation  of  statistics  -  The  SPSS  regression  package 
allows  calculation  of  numerous  regression  coefficients, 
statistics  and  residual  statistics  (difference  between 
observed  and  calculated  values).  Also  possible  are  scatter 
plots  of  residuals  and  statistical  tests  for  residual 
analysis,  etc. 

The  regression  constant  and  regression  coefficients  are  determined 
in  order  to  minimize  the  sum  of  the  squared  residuals.  The  residu- 
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als  are  the  difference  between  the  observed,  dependent  variable  and 
the  prediction  by  the  regression  equation.  The  SPSS  program  auto- 
matically  includes  those  independent  variables  which  meet  the  957. 
confidence  level  based  on  the  computed  value  of  the  F  statistic. 
Only  those  variables  which  meet  the  specified  level  at  any  given 
step  are  retained  in  the  regression  equation  and  all  those  vari- 
ables which  fall  below  the  specified  level  are  deleted  from  the 
regression  equation. 

Preliminary  regression  equations  were  developed  in  order  to  predict 
the  7Q2,  7Q20,  7050,  3Q2,  3020,  ^QsQ,  30Q2,  30020,  ^OQsO  and  O95 
low  flows  as  a  function  of  basin  physiographic  and  hydrometeorolog- 
ic  parameters.  Untransformed  parameters  which  entered  into  the 
equation  in  order  of  significance  are  drainage  area,  base  flow 
index,  slope  and  mean  annual  precipitation  or  snowfall.  The  resul- 
ting equation  for  7O20  and  O95  had  r2  values  of  .84,  .88,  and  a 
standard  error  of  0.50  and  0.67  respectively  (see  Table  4.1). 

4.2.3  Transformations 

Several  methods  of  transforming  the  parameters  to  normalize  the 
distribution  were  considered,  including: 

i)  LOGio  V)  ROOTS 

ii)  LOG  NORMAL  vi)  SQUARES  ) 

iii)  10^  vii)  CUBICS  )POWERS 

iv)  e^  viii)  COMBINATIONS 

Logarithmic  transformations  were  found  to  be  unsuccessful  in  the 
present  investigation. 

Since  the  logarithmetic  transformations  were  unsuccessful,  expon- 
ents (e^,  10^)  were  considered  next.  These  transformations  in- 
creased skewness  and  no  improvement  was  found  in  the  regression 
equation. 
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Plotting  of  each  independent  parameter  against  the  dependent  param- 
eter suggested  that  an  exponential  conversion  of  some  of  the  inde- 
pendent parameters  would  provide  a  better  relationship. 

The  literature  review  also  revealed  several  forms  of  prediction 
equations  using  the  square  root  of  selected  independent  parameters, 
and  hence  this  type  of  transformation  was  attempted.  The  results 
are  tabulated  in  Table  4.2.  The  r2  values  are  lower  and  there  are 
higher  standard  errors  of  SE  values  when  compared  to  equations 
developed  using  the  untransformed  parameters. 

Exponents  ranging  from  1/3  to  3  transforms  were  used  for  DA,  BFI, 
and  LEN.  The  best  general  improvement  was  found  for  DA^,  BFI^  and 
LEn2. 

Recent  literature  also  suggested  combinations  of  Independent  param- 
eters. For  example,  it  was  suggested  that: 


LEN  y 

fSLP)^ 


and 


BFI 


(len/(slp: 


(4.5; 


would  provide  a  combination  parameter  which  would  increase  the  fit 
of  the  regression  equation.  However,  it  was  found  that  these  com- 
bination transforms  did  not  significantly  improve  the  regression 

equations  and  could  possibly  result  in  spurious  results  when  com- 
bined in  this  fashion  with  drainage  area  and/or  Base  Flow  Index. 


Table  4.3  surmiarizes  the  final  transformed  and  untransformed  simple 
correlations  between  the  independent  variables  and  using  various 
transformations  on  the  independent  variables.  With  reference  to 
Table  4.3,  it  is  evident  that  DA^,  LENGTh2  and  BFl2  have  correla- 
tion coefficients  of  0.95,  0.88  and  0.31  compared  to  the  correla- 
tion coefficients  of  0.84,  .68  and  .31  for  their  untransformed 
counterparts. 
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Table  4.3 
Correlation  Hatrii  of  Dependent  and  Independent  Parateters 


R6 

NAP 

HAS 

HAR 

DA 

ACLS 

LAT 

L0N6 

LEN 

SLP 

RECORD 

BFI 

REG 

(n) 

(») 

(u) 

(ki'2) 

(X) 

(degree) 

(decree) 

(ki) 

(■/■) 

(years) 

(index) 

(code: 

R8 

1.00 

-O.20 

-0.48 

-0.32 

-0.37 

-0.13 

-0.32 

-0.65 

-0.38 

0.26 

-0.10 

0.31 

0.33 

NAP 

-0.20 

1.00 

0.55 

0.44 

0.11 

0.15 

0.54 

-0.20 

0.15 

-0.03 

0.19 

0.32 

0.14 

HAS 

-0.48 

0.55 

1.00 

0.57 

0.25 

0.42 

0.88 

-0.05 

0.18 

0.27 

0.18 

0.07 

-0.07 

HAR 

-0.32 

0.44 

0.57 

1.00 

0.19 

0.18 

0.48 

-0.02 

0.14 

0.06 

0.16 

0.20 

0.00 

DA 

-0.37 

0.11 

0.25 

0.19 

1.00 

-0.05 

0.21 

0.19 

0.82 

-0.28 

0.58 

0.11 

-0.14 

ACLS 

-0.13 

0.15 

0.42 

0.18 

-0.05 

1.00 

0.46 

0.03 

-0.05 

0.06 

0.04 

0.18 

-0.15 

LAT 

-0.32 

0.54 

0.88 

0.48 

0.21 

0.46 

1.00 

-0.22 

0.16 

0.27 

0.14 

0.22 

-0.07 

LONG 

-0.65 

-0.20 

-0.05 

-0.02 

0.19 

0.03 

-0.22 

1.00 

0.15 

-0.26 

-0.12 

-0.35 

-0.34 

LEN 

-0.38 

0.15 

0.18 

0.14 

0.82 

-0.05 

0.16 

0.15 

1.00 

-0.44 

0.45 

0.05 

-0.17 

SLP 

0.26 

-0.03 

0.27 

0.06 

-0.28 

0.06 

0.27 

-0.26 

-0.44 

1.00 

-0.18 

0.19 

0.06 

RECORD 

H).10 

0.19 

0.18 

0.16 

0.58 

0.04 

0.14 

-0.12 

0.45 

-0.18 

1.00 

0.19 

0.04 

BFI 

0.31 

0.32 

0.07 

0.20 

0.11 

0.18 

0.22 

-0.35 

0.05 

0.19 

0.19 

1.00 

0.20 

RES 

0.33 

0.14 

-0.07 

0.00 

-0.14 

-0.15 

-0.07 

-0.34 

-0.17 

0.06 

0.04 

0.20 

1.00 

7Q2 

-0.17 

0.11 

0.20 

0.15 

0.89 

-0.01 

0.24 

0.05 

0.70 

-0.11 

0.57 

0.32 

-0.06 

7820 

-0.16 

0.10 

0.18 

0.13 

0.84 

-0.02 

0.24 

0.05 

0.68 

-O.09 

0.50 

0.31 

-O.07 

7850 

-0.15 

0.09 

0.18 

0.12 

0.83 

-0.02 

0.24 

0.05 

0.67 

-0.09 

0.48 

0.31 

-0.07 

3Q2 

-0.17 

0.11 

0.20 

0.15 

0.89 

-0.01 

0.24 

0.05 

0.70 

-0.11 

0.57 

0.31 

-0.07 

3B20 

-0.16 

0.09 

0.18 

0.13 

0.84 

-0.02 

0.24 

0.05 

0.68 

-0.09 

0.50 

0.30 

-0.08 

3850 

-0.15 

0.09 

0.18 

0.13 

0.83 

-0.02 

0.24 

0.05 

0.67 

-0.09 

0.47 

0.30 

-O.07 

30S2 

-0.16 

0.13 

0.19 

0.15 

0.89 

-0.01 

0.23 

0.04 

0.71 

-0.13 

0.56 

0.34 

-0.04 

30880 

-0.13 

0.11 

0.16 

0.11 

0.84 

-0.03 

0.22 

0.03 

0.68 

-0.11 

0.49 

0.34 

-0.03 

30850 

-0.12 

0.10 

0.15 
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4.2.4  Regional  Regression  Results 

The  form  of  the  equation  of  the  final  regressions,  the  correspond- 
ing coefficients,  etc.  are  shown  in  Table  4.4.  (The  square  root  of 
the  period  of  record  was  used  as  a  weighting  factor  for  each  sta- 
tion.) 

It  should  be  noted  that  DA^,  BFI^,  and  LEN^  entered  sequentially 
into  the  equation.  The  length  parameter  was  not  significant  for 
30-day  average  low  flow  values  with  a  recurrence  interval  greater 
than  2  years.  The  constants  of  the  equation  are  all  relatively 
small.  This  is  due  to  the  large  values  created  by  transforming  the 
gauged  data  to  exponent  values  of  cubed  and  squared  (e.g.  mean  DA 
is  414  km2,  and  the  cubic  transform  is  7.1  x  10^).  Hence  the  small 
coefficients  are  statistically  significant  due  to  the  large  inde- 
pendent parameters  with  which  they  are  associated.  The  coeffic- 
ients decrease  with  increased  recurrence  interval  and  increase  with 
the  increase  of  the  n-day  flows  analysed. 

4.2.5  Sub-regions 

It  is  possible  that  sub-regions  could  be  identified  for  which  local 
regional  equations  may  result  in  a  better  fit  to  the  data  since  any 
extraneous  factors  affecting  the  low  flow  which  are  not  included  in 
the  regression  would  be  nearly  constant  in  the  region.  In  other 
words,  regional ization  in  this  context  is  an  attempt  to  account  for 
neglected  factors,  that  is  hydrologic  factors  which  may  be  common 
to  a  region.  A  different  set  of  prediction  equations  is  developed 
for  identified  sub-regions  which  indirectly  takes  into  account 
local  low  flow  characteristics. 

A  preliminary  regression  equation  was  developed  covering  all  sta- 
tions (i.e.  no  regions  initially,  see  Table  4.4).  An  analysis  of 
residuals  was  then  undertaken  by  plotting  the  standardized  resid- 
uals on  a  colour  coded  map  of  the  hydrometric  gauge  locations  (see 
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Figure  4.1).  From  an  analysis  of  the  plotted  residuals,  it  was 
possible  to  identify  a  cluster  of  residuals  which  indicated  that 
the  preliminary  regression  equation  consistently  overestimated  the 
low  flow  estimates  in  the  streams  draining  to  Lake  Huron  and  East 
Lake  Erie.  Similarly,  a  cluster  of  residuals  around  West  Lake  Erie 
consistently  underestimated  the  low  flows. 

Subsequent  to  the  preliminary  identification  of  low  flow  regions, 
separate  prediction  equations  were  formulated  for  each  of  the  three 
sub-regions  referred  to  above,  that  is  namely: 

1.  Lake  Huron 

2.  West  Lake  Erie 

3.  East  Lake  Erie 

Refer  to  Figure  4.1  for  identification  of  boundaries  locating  each 
of  these  regions. 

The  standard  error  of  estimate  of  the  prediction  equations  for  each 
region  are  considerably  better  than  the  standard  error  of  estimate 
for  the  prediction  equation  developed  for  the  Southwestern  West 
Central  Region  as  a  whole  (see  Tables  4.4  and  4.5  for  comparison  of 
standard  error  of  estimates  of  the  region  and  sub-regions  predic- 
tion equations). 

A  regional  set  of  regression  equations  were  developed  for  the  three 
regions  indicated  earlier  (see  Figure  4.1)  for  which  the  coeffic- 
ients and  equation  statistics  are  summarized  in  Table  4.5. 

In  the  Lake  Huron  region  only  DA^  and  BFI^  entered  into  most  of  the 
equations,  while  in  some  cases  DA^  was  the  only  significant  para- 
meter. The  West  Lake  Erie  region  had  Log  (MAR)  enter  into  5  of  the 
extreme  value  equations,  BFI^  and  LEN^  remained  as  significant 
parameters  and  DA^  did  not  appear  to  be  significant  except  for 
short  recurrence  interval  low  flows.  Alternately  Log  (MAS)  entered 
into  the  regression  equations  for  the  East  Lake  Erie  region  while 
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in  some  cases  BFI^  did  not  appear  significant.  In  the  Lake  Huron 
region  DA^  and  LEN^  become  significantly  intercorrelated  and  hence 
LEn2  is  omitted  from  the  resulting  equation.  Stream  length  LEN^ 
appears  to  explain  the  variation  of  low  flow  values  better  than  DA^ 
and  hence  due  to  significant  intercorrelation  DA^  is  omitted  from 
the  resulting  equations  for  West  Lake  Erie  region.  Base  flow  index 
BFl2  parameter  does  not  vary  areally  within  the  East  Lake  Erie 
region  and  hence  does  not  explain  the  low  flow  value  variation  in 
this  region  significantly. 

Although  R  and  SE  values  improved,  they  must  be  examined  in  respect 
to  the  diminished  samples  from  which  they  were  determined.  Con- 
stants for  equations  for  which  BFI^  or  DA-^  were  not  included  tend 
to  be  large  positive  values  when  compared  to  equations  with  these 
parameters. 

4.2.6  Summary  of  Regression  Analysis 

•  Drainage  Area  (DA),  Base  Flow  Index  (BFI)  and  Length  (LEN)  were 
found  to  be  the  main  physiographic  parameters  which  affect  low 
flows. 

•  Mean  Annual  Snowfall  (MAS)  and  Mean  Annual  Runoff  (MAR)  entered 
prediction  equations  for  sub-regions.  Figure  4.1  depicts  the 
sub-regional  boundaries  (Lake  Huron,  West  Lake  Erie  and  East 
Lake  Erie) 

•  Table  4.4  summarizes  the  recommended  regression  equations  for 
six  average  d-day  flow  for  various  recurrence  interval  low 
flows 

•  Transforming  drainage  area  by  exponent  3,  BFI  and  length  by 
exponent  2  and  by  taking  the  Log^o  of   MAR  and  MAS  produced 
better  regression  results 

•  Table  4.5  summarizes  the  regression  equations  for  the  three 
sub-regions. 
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4.3    Index  Method 

The  previous  analysis  identified  possible  sub-regions  and  determin- 
ed that  drainage  area  (DA)  is  highly  correlated  to  low  flows. 

Therefore,  a  graph  of  7Q2  as  a  function  of  DA  was  plotted  (see 
Figure  4.2).  This  relationship  had  an  r2  of  .48  and  SE  of  .39  for 
a  sample  size  of  64  stations.  Figure  4.3  illustrates  values  of  O95 
plotted  against  drainage  area. 

Three  stations;  02FC013,  North  Saugeen  River  near  Paisley;  02GC007, 
Big  Creek  near  Walsingham;  and  02GC026,  Big  Otter  Creek  near 
Colton,  exhibit  abnormally  high  low  flows.  Subsequent  investiga- 
tion identified  that  these  stations  have  significantly  higher  BFI 
values.  Hence  some  areas  require  more  than  contributing  drainage 
area  for  good  estimates  of  low  flow  characteristics  (see  Section 
4.5). 

Further  subdivided  into  sub-regions,  Figure  4.4  illustrates  how  the 
7Q2  flow  varies  per  sub-region  as  a  function  of  DA.  The  subregion- 
al  analysis  indicates  a  minimum  requirement  of  drainage  area  of  30 
to  150  km2  for  non-zero  7Q2  values  in  Lake  Huron  and  East  Lake  Erie 
regions. 

Previous  studies  have  shown  interrelationships  between  n  day  and  10 
day  average  low  flows.  Further  to  this,  the  ratio  of  various  n  day 
flows  with  y  year  recurrence  to  7Qy  was  determined.  The  results  of 
these  analyses  are  depicted  on  Figures  4.5  and  4.6. 

To  use  this  method  knowing  the  drainage  area,  the  7Q2  flow  could  be 
determined  from  Figure  4.2  or  Figure  4.4.  If  another  nQy  flow  is 
required,  it  can  be  estimated  by  using  Figure  4.5  and/or  Figure  4.6 
to  determine  the  appropriate  flow  ratio  as  a  function  of  7Q2.  The 
estimating  error  for  7Q2  and  the  errors  associated  with  the  index 
ratio  suggest  that  this  method  would  not  provide  the  same  level  of 
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7Q2  AS  A  FUNCTION  OF  DRAINAGE  AREA 
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accuracy  as  the  regression  equations.  This  is  examined  in  more 
detail  in  Section  5.0. 

The  graphs  of  n-day  low  flows  to  7-day  low  flows  as  well  as  the 
graph  of  7-day  y  year  recurrence  ratio  to  7Q2  low  flows  give  some 
insight  to  the  interrelationships  between  the  various  extreme  value 
low  flows.  For  example,  the  ratio  of  30  day  average  low  flows  to  7 
day  low  flows  tends  to  increase  with  increased  recurrence  interval 
while  the  ratio  of  3  day  average  low  flows  to  7  day  average  low 
flows  decrease  with  increased  recurrence  interval.  As  expected, 
the  30  day  low  flows  are  greater  than  7  day  low  flows  (i.e.  ratio 
ranges  from  1.5  to  1.8  times  the  7  day  low  flow  values).  Figure 
4.5  illustrates  the  typical  ratios  of  7  day  y  recurrence  low  flows 
as  a  ratio  of  7  day  2  year  recurrence  low  flow  and  how  the  ratio 
decreases  with  increased  recurrence  interval. 

It  may  be  possible  to  relate  shorter  periods  of  record  to  7Q2  flows 
and  then  use  the  relationships  developed  here  to  estimate  longer 
recurrence  low  flows.  This  concept  is  recommended  for  further 
research,  since  it  is  beyond  the  scope  of  this  investigation. 

4.4    Mapping  Proration 

In  the  past  it  has  been  common  practice  to  prorate  unit  flows  from 
nearby  gauged  watersheds  to  establish  the  required  low  flows.  This 
is  done  by  experienced  personnel  who  have  a  good  understanding  of 
local  stream  characteristics  and  other  factors  within  the  region 
(i.e.  diversions  and  regulation,  etc.). 

For  this  analysis,  the  unit  area  7Q20  and  Q95  flows  (l/s/km2)  were 
summarized  (as  per  Section  3.2.7)  for  stations  in  the  region  (see 
Figure  4.7,  in  pocket). 

The  methods  of  reciprocal  distance  percentage  could  be  used  in  the 
proration  estimation  (i.e.  the  distance  to  nearby  gauges  would  be 
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estimated  and  unit  values  could  then  be  weighted  by  the  reciprocal 
of  the  distance  as  a  percentage  and  then  averaged). 

A  statistical  summary  of  average  unit  low  flow  characteristics  was 
also  undertaken  and  is  summarized  in  Table  4.6.  Area  average  unit 
low  flow  rates  (l/s/km^)  together  with  standard  deviation  for  the  5 
durations  were  tabulated  for  the  region  for  the  2  to  50  year  recur- 
rence intervals.  It  is  interesting  to  note  that  in  most  cases,  the 
standard  deviations  are  approximately  equal  to  or  greater  than  the 
corresponding  mean  flows. 

4.5    Mapping  Low  Flows 

From  Section  4.2  it  was  determined  that  BFI  was  the  second  most  im- 
portant parameter  in  predicting  low  flows.  Moin  and  Shaw,  1986 
have  recently  mapped  isolines  of  BFI  on  a  map  of  Southern  Ontario. 
If  it  is  theoretically  possible  to  map  BFI  (assumed  to  be  a  measure 
of  basin  geology  and  base  flow  characteristics),  then  it  should  be 
possible  to  map  unit  area  low  flows  due  to  the  high  correlation 
previously  noted.  A  regression  equaton  relating  BFI  only  to  unit 
7020  flows  was  determined.  The  equation  was  applied  to  mapped 
isolines  of  BFI  and  the  BFI  isolines  were  transferred  into  unit 
7Q20  isolines  (see  Figure  4.8).  Hence,  by  identifying  the  site 
location  on  the  map  estimating  the  unit  7Q20  from  the  isolines  and 
multiplying  by  the  drainage  area,  7Q20  estimates  for  ungauged  sites 
could  be  determined. 

The  r2  was  0.68  and  0.53  with  S.E.  of  0.89  and  0.75  for  regression 
equations  of  7Q20  and  Qgs  as  a  function  of  BFI.  (For  this  analy- 
sis, 7O20  and  Q95  were  converted  to  unit  values  per  kilometer 
squared  and  from  cubic  metres  to  litres.)  The  resulting  equation 
illustrated  that  areas  with  BFI  less  than  0.2  would  have  7020 
values  of  zero. 
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The  method  is  tested  as  described  in  Section  5.0.  Further  investi- 
gation should  examine  the  contouring  of  unit  low  flow  values 
directly  for  comparison  of  this  technique. 

4.6    Summary  of  Estimation  Techniques 

The  regression  equations  summarized  in  Table  4.2  appear  to  best 
represent  low  flows  (S.E.  .19  to  .41).  Other  studies  (see  Section 
2.0)  found  the  S.E.  of  prediction  to  range  from  .30  to  .37  for 
similar  sample  sizes  (Condie,  1985,  Chang,  1977,  Logan,  1986). 
Therefore,  the  prediction  errors  determined  in  the  present  investi- 
gation are  of  smaller  magnitude. 

The  standard  error  for  the  BFI  isoline  conversion  equation  is  high- 
er due  to  additional  error  introduced  by  transforming  7Q20  and  Qgs 
by  division  of  drainage  area  and  conversion  from  m^/s  to  1/s. 
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5.0  TESTING  REGIONALIZATION  METHOD 

5.1  General 

The  relative  prediction  accuracy  of  the  four  methods  developed  in 
Section  4.0  was  examined  using  10  stations  not  included  in  the 
development  of  the  regional  low  flow  estimation  techniques. 

The  10  test  stations  all  have  at  least  10  years  of  record,  and  are 
active  within  the  last  5  years.  Test  stations  which  had  indepen- 
dent parameters  outside  those  use  in  the  analysis  were  not  used 
(i.e.  02HA003,  Niagara  River  at  Queenston  [DA  =  686,000  km2]  and 
02GA032  O.A.C.  Farm  Gauge  No.  5  at  Guelph  [DA  =  2.5  km2]).  Unfor- 
tunately W.S.C.  station  02HB015,  Spencer  Creek  near  Westover  was 
subsequently  determined  to  be  heavily  regulated.  Although  similar 
to  the  analysed  data  base,  the  range  of  drainage  areas  tend  to  be 
on  the  low  side  and,  therefore,  further  testing  should  be  completed 
using  more  stations  with  a  greater  range  of  independent  parameters. 

A  list  of  the  test  stations  and  their  appropriate  independent  and 
dependent  parameters  are  given  in  Table  5.1. 

5.2  Testing 

Figures  5.1  and  5.2  summarize  the  results  of  the  analysis  tech- 
niques compared  to  estimates  from  statistically  fitted  gauged  low 
flows.  The  Nash-Sutcl iffe  1970  model  efficiency,  r2,  was  used  to 
help  assess  the  performance  of  the  identified  estimation  methods, 
as  calculated  below: 


r2  =  1  _  S(Os  -  Qo)^ 
Z(Qo  -  Qm)^ 


where  Oo  and  Qs  are  the  observed  and  simulated  discharges,  and  Qm 
is  the  mean  of  the  observed  discharges. 
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The  results  are  tabulated  in  Appendix  B.  Due  to  the  small  sample 
of  test  stations  the  comparison  results  are  used  here  only  to  iden- 
tify the  ranking  of  low  flow  estimation  of  methods  and  not  to  infer 
any  statistical  significance. 

The  "whole  area"  regression  equation  (not  sub-regions)  provided  the 
best  comparison  of  single  station  estimate  versus  "regional  method" 
estimate.  However,  this  method  tended  to  overestimate  the  majority 
of  the  7Q20  and  O95  low  flows.  Further  analysis  (i.e.  use  of  sub- 
regional  equation)  did  not  give  comparable  results  when  applied  to 
the  test  stations. 

The  remaining  three  methods  ranked  2,  3,  4  as  Isoline  Method,  Index 
Method,  Proration  Method  respectively. 

The  isoline  method  resulted  in  fairly  good  results  but  should  be 
tested  with  a  larger  range  and  number  of  stations.  This  method 
tended  to  overestimate  smaller  low  flows  and  underestimate  larger 
low  flows  in  this  sample  test  data  set.  The  heavily  regulated  sta- 
tion was  estimated  by  the  isoline  method  better  than  the  regression 
equations  which  were  formed  without  heavily  regulated  stations. 
Further  investigation  into  mapping  isolines  of  low  flow  should, 
therefore,  be  undertaken.  This  method  had  the  highest  Nash- 
Sutcliffe  r2  value  for  estimating  Q95  for  the  10  test  stations. 

The  index  method  showed  similar  results  to  that  of  the  isolines 
method  in  that  it  overestimated  smaller  values  and  underestimated 

larger  values. 

The  proration  method  (using  nearby  gauges)  did  not  perform  as  well 
as  the  three  alternate  estimation  methodologies  developed  from  this 
study. 

Therefore,  as  an  overall  conclusion,  the  Regression  technique  (hav- 
ing a  Nash  Sutcliffe  r2  of  .69)  relating  drainage  area,  base  flow 
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index  and  length  to  low  flows  has  been  determined  as  the  best 
method  for  estimating  7Q20  low  flows  for  non-regulated  ungauged 
watersheds  with  drainage  areas  14  km^  to  4000  km^,  BFI  indices  .10 
to  .80  and  Lengths  from  6  km  to  190  km. 

5.3    Sensitivity  Analysis 

The  Regression  equations  identified  above  for  7O20  and  Q95  values 
were  examined  with  respect  to  sensitivity  of  parameter  estimation. 
Table  5.2  suirmarizes  the  percentage  difference  from  the  mean  para- 
meter estimation  by  varying  each  parameter  by  ±15%.  It  is  noted 
that  7Q20  and  Q95  values  are  sensitive  to  BFI  and  Length  parameters 
(15%  change  causes  27  to  38%  and  15  to  25%  change  for  BFI  and  LEN 
respectively  in  the  estimated  value).  Both  BFI  and  LEN  overesti- 
mates cause  higher  errors  than  the  corresponding  underestimates. 
Therefore,  care  should  be  taken  in  estimating  these  parameters  for 
ungauged  watersheds. 


TABLE  5.2 
Surmiary  of  Sensitivity  Analysis 

7Q20  =  -.166  +  9.03x10-11  qa^  +  1.10  BFl2  +  8.35x10-5  len2 
O95  =  -.23  +  1.25x10-10  da3  +  1.98  BFI^  +  8.97x10-5  LEn2 

Values  of  7O20  (m^/s)  Varying  with  ±15%  in  Indicated  Parameters 


DA 

BFI 

LEN 

7Q20 
(m3/s) 

(%  Diff) 

7Q2O 
(m3/s) 

(%  Diff) 

7Q20 
(m^/s) 

(%  Diff) 

.227 
.230 
.233 

7 
0 
7 

.167 
.230 
.303 

27 

0 

30 

.180 
.230 
.288 

21 

0 

25 

Values  of  O95  (m^/s)  with  Varying  Parameters 


DA 

BFI 

LEN 

995 
(m3/s) 

(%  Diff) 

Q95 
(m^/s) 

(%  Diff) 

Q95 
(m3/s) 

(%  Diff) 

.348 

1 

.240 

31 

.299 

15 

.352 

0 

.352 

0 

.352 

0 

.356 

1 

.484 

38 

.414 

18 
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6.0  CONCLUSIONS  AND  RECOMMENDATIONS 

6.1  Conclusions 

1.  Regressions  utilizing  Drainage  Area,  Base  Flow  Index  and 
Stream  Length  as  independent  parameters  provide  reasonable 
estimates  of  low  flow  values  (see  Table  4.4). 

2.  The  equations  for  estimating  low  flows  are  sensitive  to  vari- 
ations in  values  of  Base  Flow  Index  (BFI)  and  Length  (LEN). 
Care  should  be  taken  in  estimating  these  parameters  for 
ungauged  watersheds.  From  the  literature  review,  it  is 
apparent  that  even  short  term  gauging  (i.e.  2  years)  can  pro- 
vide good  estimates  of  BFI.  Alternatively,  available  maps  of 
BFI  isolines  can  be  used. 

3.  The  index,  or  graphical  method,  was  also  found  to  give 
reasonable  estimates  of  low  flow  characteristics.  Both  the 
regression  and  graphical  methods  appear  to  provide  low  flow 
estimates  superior  to  those  obtained  solely  from  area  pro- 
ration to  gauged  sites. 

4.  Statistical  tests  on  the  data  base  indicated  that  the  series 
of  extreme  low  flows  may  exhibit  some  trend  and  persistence 
with  some  possibility  of  non-random  characteristics. 

5.2  Recommendations 

1.  The  reasons  for  trend,  persistence  and  non-randomness  charac- 
teristics should  be  further  examined.  This  should  include  an 
analysis  of  the  applicability  of  the  statistical  tests  in 
addition  to  an  evaluation  of  other  factors  such  as  climatic 
trends  and  variations  in  groundwater  storage  and  aquifer  con- 
ditions. 
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2.  To  maintain  a  relevant  data  base  of  recent  low  flow  charac- 
teristics for  the  regions,  the  analysis  of  single  station  low 
flow  characteristics  should  be  updated  and  results  remapped 
every  5  years. 

3.  The  regional ization  methods  which  were  developed  should  be 
tested  on  a  larger  data  base  to  support  the  findings  of  this 
study  and  provide  refinements  to  estimating  procedures  for 
other  regions.  This  should  include  further  investigation  on 
the  form  of  the  regression  equation  and  identification  of 
sub-regions,  etc.  The  practical  advantages  of  utilizing 
untransformed  independent  parameters  (DA,  BFI,  LEN)  should  be 
further  assessed  in  order  of  scatter  at  test  locations  and 
the  small  size  of  the  regression  coefficients,  etc. 

4.  The  index  method  should  be  extended  using  all  available  short 
term  station  data  to  refine  the  estimates  of  7Q2  values  as  a 
function  of  drainage  area.  This  would  maximize  use  of  the 
available  data  base. 

5.  Isoline  maps  of  unit  low  flow  values  should  be  prepared  using 
station  estimates  and  compared  to  the  results  of  mapping 
using  the  BFI  regression/conversion  method. 

5.  These  techniques  for  estimating  low  flow  values  at  ungauged 
sites  should  be  extended  to  the  remaining  four  Ministry  of 
Environment  Regions: 

i)  Central  Ontario  Region 

ii)  Southeastern  Region 

iii)  Northwestern  Region 

iv)  Northeastern  Region 
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APPENDIX  A.l 
EXTREME  VALUE  ANALYSIS 


Introduction 

Statistical  distributions  allow  estimates  of  probability  of  exceedance  of 
events  to  be  made  by  analytical  techniques.  Various  methods  are  avail- 
able, and  a  good  discussion  is  given  by  McMahon  and  Mein  (ref  9). 

At  a  gauged  location,  average  low  flows  can  be  determined  for  selected 
durations  of  time  for  each  year  of  record.  For  example,  the  corresponding 
minimum  average  consecutive  15-day  low  flow  can  be  determined  for  each 
annual  record.  The  series  of  15-day  average  annual  low  flows  can  then  be 
ranked  lowest  to  highest  and  an  extreme  value  analysis  undertaken  using  a 
theoretical  distribution. 

The  Gumbel  Extreme  Value  Distribution  is  commonly  used  for  fitting  low 
flow  frequency  curves  to  the  available  data. 

The  Gumbel  Type  III  Distribution 

This  is  a  variation  of  Fisher  and  Tippett's  third  asymnptotic  distribution 
of  extreme  values  and  is  sometimes  referred  to  as  the  Weibull  distribu- 
tion. In  view  of  his  many  contributions  it  is  often  referrred  to  in 
hydrology  as  the  Gumbell  II.  The  probability  density  function  is: 


0{x) 


a-1 


exp 


X  -  e 

u  -  e 


(1) 


where  e  is  the  lower  bound  parameter 
u  is  the  characteristic  drought 
and   a  is  the  shape  parameter 


Various  methods  are  available  for  determining  the  distribution  parameters 
for  a  particular  data  set.  See  references  2  and  6  for  example. 
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The  Probability  Function 


The  density  function  is  integrable  and  gives  the  distribution  function 


F  (x)  =  1  -  exp 


X  -  e 

u  -  S 


(2) 


which  gives  the  probability  of  non-exceedance  of  x. 

Since  it  is  more  common  to  require  x  for  a  given  probability  of 

non-exceedance,  a  simple  re-arrangement  gives 


x  =  e 


+  (u  -  e)  j  -  In  Pi  -  F(x)  ~| 


1/a 


(3) 


There  are  occasions  when  the  sample  series  of  low  flows  can  have  a  large 
negative  coefficient  of  skewness,  and  since  the  Gumbel  Type  III  distribu- 
tion cannot  have  a  skewness  of  less  than  -1.14,  then  the  distribution 
cannot  be  fitted. 

There  are  insuffucient  natural  samples  available  with  such  low  skewness  to 
make  a  firm  choice  of  an  alternative  treatment,  but  from  the  few  avail- 
able, it  was  found  that  the  negatively  skewed  three-parameter  lognormal 
distribution  provided  an  acceptable  alternative. 


If  y  =  In  (a-x)  is  normally  distributed,  then  the  probability  density 
function  of  x  is  given  by 


$(x) 


1 


exp 


(a-x)o  yzn 


I  "i^  [i"  (^->^)  -n  '  I 


(4) 


where  m  and  o  are  respectively  the  mean  and  standard  deviation  of  the 
series  In  (a-x)  and  a  is  an  lower  boundary  parameter.  This  form  of  the 
distribution  can  only  have  negative  skewness.  The  equation  in  the  form 
In  (x-a)  is  for  distributions  with  positive  skewness. 
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Taking  moments,  re-arranging  and  replacing  them  by  their  sample  estimates 
gives: 

k3  +  3k  -  gi  =  0  (35a)   (5) 

and  after  solving  for  k 

a  =  X  -  5/k  (37a)   (6) 

m  =  ln(-s/k)  -  l/21n(k2+l)  (38a)   (7) 

o2  =  in(k2+l)  (39a)  (8) 

Thus,  the  distribution  is  completely  defined  and  the  T-year  low  flows  can 
be  computed  from: 

Qx  =  a  -  exp(m  ot)  (9) 

Additional  details  are  discussed  in  the  user's  manual. 


Source:  Condie,  R.,  L.  C.  Cheng 

"Low  Flow  Frequency  Analysis  -  Program  LOFLOW", 
Water  Resources  Branch,  Inland  Waters  Directorate 
Environment  Canada,  Ottawa,  1987 
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APPENDIX  A. 2 
FLOW  DURATION  ANALYSIS 


A  flow  duration  curve  is  a  plot  of  the  flow  of  a  stream  against  the  per- 
cent of  time  the  indicated  flow  was  equalled  or  exceeded  during  the  period 
covered  by  the  available  flow  data.  This  curve  is  extremely  useful  for 
hydro  power  studies  and  for  characterizing  the  local  streamflow  regime. 

An  empirical  procedure  was  used  to  analyse  each  sample  discharge  record. 
Flow  duration  curves  are  derived  by  rearranging  the  available  daily  or 
monthly  flow  data  in  order  of  magnitude.  The  total  time  period  represents 
100%  of  the  time.  Therefore,  by  definition  of  the  procedure  used  here, 
the  largest  value  is  exceeded  0  percent  of  the  time  and  the  smallest  value 
is  exceeded  100%  of  the  time. 

The  flow  duration  curve  represents  the  percent  of  time  that  a  specific 
discharge  occurs  at  that  location.  However,  the  curve  does  not  indicate 
the  period  of  time  in  the  year  when  the  flow  is  less  than  or  equal  to  the 
selected  value.  Therefore,  in  some  instances,  it  is  also  useful  to  devel- 
op flow  duration  curves  on  a  monthly  basis,  that  is,  all  data  for  the 
month  of  January  over  the  entire  period  of  record  is  analysed  independent- 
ly to  produce  a  flow  duration  curve  representative  of  flow  conditions  in 
January. 

Single  station  flow  duration  analyses  and  corresponding  monthly  flow  dura- 
tion curves  have  been  determined  at  all  relevant  stream  and  flow  locations 

across  the  Province. 
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APPENDIX  A. 3 

NONPARAMETRIC  TESTS  FOR  INDEPENDENCE, 

TREND  AND  RANDOMNESS 


A. 3.1    Test  Description 


This  appendix  briefly  summarizes  the  functions  evaluated  in  the 
LFA  package  and  gives  the  methods  used  to  determine  their 
statistical  significance. 


A. 3. 1.1  Introduction 


Any  statistical  test  of  significance  will  generally  be  made  using 
the  following  steps: 

a)  State  the  null  hyposthesis,  Hq.  For  instance  in  split  sample 
tests,  the  null  hypothesis  may  be  that  there  is  no  difference 
between  the  sample  means. 

b)  Choose  a  significance  level,  a. 

c)  Choose  an  appropriate  statistical  test.  In  this  program  all 
tests  are  nonparametric. 

d)  Compute  the  test  statistic. 

e)  The  sampling  distribution  of  the  test  statistic  is  known  and 
has  been  tabulated,  and  the  chosen  significance  level  then 
defines  the  region  of  rejection. 

f)  If  the  computed  test  statistic  lies  in  the  region  of 
rejection,  then  the  null  hypothesis  is  rejected. 

Consider  now  the  four  tests  in  this  program. 
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A. 3. 1.2  The  Spearman  Rank  Order  Serial  Correlation  Coefficient  for 
Independence 

If  the  series  Q.  with  i  ranging  from  1  to  N  is  put  in 
chronological  order,  two  time  series  are  formed  and  their 
respective  ranks  computed: 

Ol,  02  Qn-1  by  X.,  the  rank  of  Q. ;  i  =  l;  N-1 

and  02,  O3  On  by  y . ,  the  rank  of  0.;  1=2,  N 

then  Spearman  rank  order  serial  correlation  coefficient  is 

S^=  I    (Zx2+  2y2  -  Sdf)  (Sxf  2y2)-5^  (1) 

where  Sx^  =  (m3-m)/12-  ZTx 

2y^  =  (m3-m)/12-  ZT 

d.  is  the  difference  in  rank  between  x.  and  y. 
m  =  N-1 


and  the  summations  are  over  the  m  pairs  of  x.,  y.. 


Ignoring  for  the  moment  the  terms  in  T  and  putting  them  at  zero, 
equation  (a)  becomes 

S^=  l-(6Zd2)  (m3-m)  (2) 


the  more  familiar  form  of  the  Spearman  rank  correlation 

coefficient. 

The  terms  in  T  adjust  for  tied  ranks  and  are  computed  as 
follows.  If  for  instance  three  observations  in  the  x  series  were 
tied  for  ranks  17,  18,  and  19  then  each  observation  is  given  the 
rank  18;  if  two  were  tied  for  ranks  24  and  25,  then  each  is 
ranked  24.5. 
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For  each  tied  set,  T  is  computed  from 


T^  =  (j3  -  J)/12 


where  J  is  the  number  of  observations  tied  at  a  give  rank. 

ST  and  2T  are  defined  by  the  extension  of  the  foregoing. 
X      y 

When  N  is  10  or  greater,  then  the  function 

t  =  S^  [(m-2)/(l-S2)]^  (3) 

is  distributed  like  Student's  t  with  m-2  degrees  of  freedom.  A 
one-tail  test  must  be  used. 

A. 3. 1.3  The  Spearman  Rank  Order  Correlation  Coefficient  Test  for  Trend 


If  the  series  Q.  with  i  ranging  from  1  to  N  is  put  in 
chronological  order,  ranks  assigned  and  denoting  the  series 


Ol,  02  Qn  by  y^,  the  rank  of  Q. 

and  1,  2,  N  by  x.,  the  sequential  order  of  Q. 


then  the  Spearman  rank  order  correlation  coefficient  r  is 
calculated  as  in  equation  1,  except  that  m  =  N,  T  =  0,  and  the 

summations  are  taken  over  the  N  pairs  of  x. ,  y.. 


For  N  =  10  or  greater,  then  the  function 

t  =  r^  [(N-2)/(l-r2)]^  (4) 

is  distributed  like  Student's  t  with  N-2  degrees  of  freedom.  The 
null  hypothesis  is  that  there  is  no  trend,  either  upward  or 
downward  with  time,  and  so  a  two-tail  test  is  used. 
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A. 3. 1.4  Runs  Above  and  Below  the  Median  for  General  Randomness 

This  randomness  test  is  based  on  the  order  or  sequence  in  which 
the  individual  scores  or  observations  were  obtained.  Actually, 
the  test  is  based  on  the  number  of  runs  which  a  sample  exhibits. 

A  run  is  defined  as  a  succession  of  identical  symbols  which  are 
followed  and  preceded  by  different  symbols  or  by  no  symbols  at 

all. 

The  total  number  of  runs  in  a  sample  of  any  given  size  gives  an 
indication  of  whether  or  not  the  sample  is  random.  If  very  few 
runs  occur,  a  time  trend  or  some  bunching  due  to  lack  of  indepen- 
dence is  suggested.  If  a  great  many  runs  occur,  systematic 
short-period  cyclical  fluctuations  seem  to  be  influencing  the 
sample. 

For  example,  once  the  median  of  the  sample  has  been  determined, 
each  observation  can  be  labelled  as  being  above  and  equal  to  or 
below  and  equal  to  the  median.  If  "A"  represents  above  and  equal 
to  the  median  and  "B"  represents  below  and  equal  to  the  median, 
then  a  sample  may  be  ordered  as 

AABBBABBBBAABA 

(A  run  represents  a  succession  of  identical  symbols).  For  our 
example,  each  run  is  underscored  and  numbered  consecutively: 


AA 

BBB 

A 

BBBB 

AA 

B 

A 

1 

2 

3 

4 

5 

6 

7 

This  sample  begins  with  2  observations  above  or  equal  to  the 
median,  followed  by  a  run  of  3  observations  below  or  equal  to  the 
median,  etc. 

Seven  runs  are  observed  in  all:  that  is,  the  total  number  of 
runs  above  and  below  the  median  RUNAB,  is  7.  If  n]  represents 
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the  number  of  events  of  type  A,  then  nj  =  5.  If  n2  denotes  the 
number  below  the  median,  type  B,  then  n2  =  8.  Thus,  the  number 
of  observations  is  equal  to  (n^  +  n2). 

In  order  to  apply  this  run  test,  one  must  determine  n^,  n2  and 
RUNAB. 


The  null  hypothesis,  H  ,  is  that  the  A's  and  B's  occur  in  random 
order.  The  alternate  hypothesis,  H. ,  is  that  the  order  of  the 
A's  and  B's  deviates  from  randomness. 


When  either  n^  or  n2  is  greater  than  20,  the  sampling  distribu- 
tion of  RUNAB  tends  to  normality  with 


RUNAB  -  [(2nin2)/(ni+n2)+l]  | 

(5) 


2        H, 

{2nin2(2nin2-ni-n2)/[(ni+n2)  (ni+n2-l)]^ 

where  z  is  an  N(0,1)  variate  as  described  in  Table  A. 4.  This 
package  uses  a  region  of  rejection  defined  by 

z  greater  than  1.96  for  a  =  .05 

z  greater  than  2.575  for  a  =  .01 

A. 3. 2    Summary  of  Tests  Results 

All  stations  for  each  region  are  tabulated  as  to  whether  the 
station  data  shows  independence,  free  from  trend  and  random  for 
significant  levels  of  1%  and  5%  in  Table  3.2,  Section  3.0).   It 
can  be  seen  from  the  summary  table  that  a  large  number  of  the 
data  sets  for  the  stations  fail  the  non-parametric  test  to  the 
point  that  the  binomial  probability  of  this  number  of  failures  is 
"0".  Therefore,  the  data  is  dependent,  is  not  free  from  trend 
and  is  not  random.  Hence,  some  degree  of  scepticism  should  be 
used  when  using  the  extreme  value  predicted  flows  from  this 
study. 
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To  examine  effects  of  length  of  record  and  regulation,  the 
analysis  was  completed  using  criteria  of  >20  years  of  record,  <20 
years  of  record  and  regulated  vs.  non-regulated  station  data. 
These  subgroups  are  also  indicated  on  Table  3.2. 

The  regulated  stations  show  a  greater  percentage  of  values  which 
can  be  said  to  pass  the  non-parametric  test  compared  to  those 
which  failed,  although  the  binomial  probability  is  still  "0". 

The  data  sets  for  stations  with  less  than  20  years  of  record 
suggest  that  neither  the  length  of  record  or  the  regulation  can 
account  for  the  dependence,  trend  and  non-randomness  of  the  data 
sets. 

It  is  also  noted  that  the  data  tended  to  be  more  significant  with 
respect  to  the  increasing  length  of  duration  for  which  the  sample 
represents.  Further  investigation  is  beyond  the  scope  of  this 
study. 


APPENDIX  B 
TEST  STATION  ANALYSIS 
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APPENDIX  B 


Table  B.l  summarizes  the  Nash-Sutcl iffe  model  efficiency  r2. 


7(0  -  Q  )^ 

2(Qo  -  V 


where  Q„  =  observed  flow 

0 

Q  =  simulated  flow 

Q  =  mean  of  observed  flow 
m 


The  efficiency  equation  compares  the  simulated  and  observed  flow 
differences  compared  to  the  observed  differences. 

The  first  part  of  Table  B.l  illustrates  the  observed  flows  and  the 
denominator  of  the  efficiency  equation.  The  following  sections 
summarize  the  resulting  estimates  via  the  four  estimation  techniques 
and  the  numerator  of  the  efficiency  equation.  Note  the  Nash- 
Sutcliffe  efficiency  equation  r2  is  not  equivalent  to  regression 
equations  coefficient  of  determination. 


TABLE  B . 1 
COMPARISON  OF  ESTIMATED  VS  ACTUAL  LOW  FLOWS   (9  STATIONS) 

STATISTICALLY  COMPUTED  GAUGED  INFORMATION 


7Q2 

7Q20 

(Q0-QM)''2 

02FD001 

0.003 

0.000 

0.004 

02GD008 

0.090 

0.007 

0.003 

02GD010 

0.030 

0.001 

0.004 

02GD019 

0.021 

0.000 

0.004 

02GE007 

0.048 

0.018 

0.002 

02HB010 

0.103 

0.034 

0.001 

02HB011 

0.455 

0.328 

0.070 

02HB012 

0.066 

0.041 

0.001 

02HB013 

0.211 

0.144 

0.006 

AVG 

0.11 

0.06 

0.01 

STD 

0.13 

0.10 

0.02 

SUM 

1.03 

0.57 

0.10 

095 

(QO-QM) 

0.000 

0.017 

0.085 

0.002 

0.028 

0.010 

0.029 

0.010 

0.065 

0.004 

0.119 

0.000 

0.531 

0.162 

0.076 

0.003 

0.227 

0.010 

0.13 

0.02 

0.16 

0.05 

1.16 

0.22 

ALL  STATIONS  REGRESSION 


7Q2 

7Q20 

(QS-QO) 

02FD001 

0.000 

0.000 

0.000 

02GD008 

0.126 

0.027 

0.000 

02GD010 

0.081 

0.002 

0.000 

02GD019 

0.000 

0.000 

0.000 

02GE007 

0.000 

0.000 

0.000 

02HB010 

0.353 

0.179 

0.021 

02HB011 

0.571 

0.324 

0.000 

02HB012 

0.165 

0.062 

0.000 

02HB013 

0.421 

0.232 

0.008 

AVG 

0.19 

0.09 

0.00 

STD 

0.20 

0.12 

0.01 

SUM 
R2 

1.72 

0.83 

0.03 
0.69 

REGIONAL  REGRESSION  EQUATIONS 


7Q2 

7Q20 

(QS-QO) "2 

02FD001 

0.000 

0.000 

0.000 

02GD008 

0.160 

0.062 

0.003 

02GD010 

0.113 

0.028 

0.001 

02GD019 

0.000 

0.000 

0.000 

02GE007 

0.000 

0.013 

0.000 

02HB010 

0.279 

0.226 

0.037 

02HB011 

0.496 

0.236 

0.009 

02HB012 

0.173 

0.221 

0.032 

02HB013 

0.399 

0.000 

0.021 

AVG 

0.18 

0.09 

0.01 

STD 

0.17 

0.10 

0.01 

SUM 

1.62 

0.78 

0.10 

R2 

-0.07 

Q95 

(QS-QO) ~2 

0.000 

0.000 

0.132 

0.002 

0.084 

0.003 

0.000 

0.001 

0.000 

0.004 

0.404 

0.081 

0.663 

0.018 

0.189 

0.013 

0.494 

0.071 

0.22 

0.02 

0.23 

0.03 

1.97 

0.19 

0.11 

Q95 

(QS-QO) ^2 

0.000 

0.000 

0.192 

0.011 

0.136 

0.012 

0.000 

0.001 

0.000 

0.004 

0.517 

0.159 

0.535 

0.000 

0.509 

0.187 

0.000 

0.052 

0.21 

0.05 

0.23 

0.07 

1.89 

0.43 

-0.96 

TABLE  B.l 
COMPARISON  OF  ESTIMATED  VS  ACTUAL  LOW  FLOWS 

INDEX  METHOD 


(9  STATIONS) 


7Q2 

7Q20 

(QS-Q0)"2 

02FD001 

0.219 

0.088 

0.008 

02GD008 

0.268 

0.107 

0.010 

02GD010 

0.215 

0.086 

0.007 

02GD019 

0.094 

0.038 

0.001 

02GE007 

0.270 

0.108 

0.008 

02HB010 

0.232 

0.093 

0.003 

02HB011 

0.305 

0.122 

0.042 

02HB012 

0.144 

0.057 

0.000 

02HB013 

0.122 

0.049 

0.009 

AVG 

0.21 

0.08 

0.01 

STD 

0.07 

0.03 

0.01 

SUM 
R2 

1.87 

0.75 

0.09 
0.06 

PRORATED  UNIT  GAUGED  VALUES 
7Q2      7Q20     (QS-Q0)~2 


I02FD001 

0. 

,236 

0.122 

0.015 

02GD008 

0. 

,236 

0.102 

0.009 

02GD010 

0. 

,122 

0.039 

0.001 

02GD019 

0. 

,058 

0.025 

0.001 

02GE007 

0. 

,242 

0.119 

0.010 

32HB010 

0. 

,131 

0.067 

0.001 

32HB011 

0. 

,196 

0.114 

0.046 

32HB012 

0. 

,055 

0.017 

0.001 

D2HB013 

0. 

032 

0.007 

0.019 

WG 

c 

).15 

0.07 

0.01 

STD 

c 

).08 

0.04 

0.01 

3UM 

] 

-.31 

0.61 

0.10 

^2 

-0.08 

ISOLINE  MAPS  FROM  BF 

7Q20 

(QS-Q0)^2 

D2FD001 

0.165 

0.027 

D2GD008 

0.070 

0.004 

)2GD010 

0.072 

0.005 

D2GD019 

0.039 

0.002 

32GE007 

0.073 

0.003 

)2HB010 

0.033 

0.000 

)2HB011 

0.209 

0.014 

)2HB012 

0.016 

0.001 

)2HB013 

0.081 

0.004 

iVG 

0.08 

0.01 

3TD 

0.06 

0.01 

5UM 

0.76 

0.06 

12 

0.37 

095 

(QS-Q0)*2 

0.312 

0.098 

0.260 

0.031 

0.140 

0.013 

0.314 

0.081 

0.277 

0.045 

0.349 

0.053 

0.189 

0.117 

0,167 

0.008 

0.169 

0.003 

0.24 

0.05 

0.07 

0.04 

2.18 

0.45 

-1.07 

Q95 

(QS-Q0)-2 

0.221 

0.049 

0.179 

0.009 

0.088 

0.004 

0.056 

0.001 

0.203 

0.019 

0.140 

0.000 

0.242 

0.084 

0.056 

0.000 

0.027 

0.040 

0.13 

0.02 

0.08 

0.03 

1.21 

0.21 

> 

0.05 

095 

(QS-Q0)~2 

0.002 

0.000 

0.177 

0.008 

0.133 

0.011 

0.012 

0.000 

0.000 

0.004 

0.372 

0.064 

0.751 

0.048 

0.139 

0.004 

0.179 

0.002 

0.20 

0.02 

0.23 

0.02 

1.76 

0.14 

0.34 

^N 
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UNIT  LOW  FLOW  CHARACTERISTICS 
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